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SUMMARY 


A method is presented for the approximate calculation of the 
stresses and distortions in a box beam representing the main structural 
component of a swept wing with a carry-through section and with bulk- 
heads perpendicular to the spars. The outer and carry-through sections 
of such a wing can be analyzed by existing methods if some means is pro- 
vided for analyzing the triangular section which connects them. In the 
method presented the triangular section is divided into free bodies and 
then equilibrium and continuity are established between these bodies. 

The result is a system of linear equations which can be solved for the 
rotations and translations of the three vertical edges of the triangular 
section. 

The application of the method is illustrated by a numerical example 
and the results are compared with previously published test data. The 
agreement is fair, with the principal discrepancies being due to the 
fact that the method is based on a very simple type of idealized struc- 
ture which prevents the appearance of shear lag in the results. Exten- 
sion of the basic approach to permit the inclusion of shear lag is 
indicated. The effects of the shear and bending flexibility of the 
bulkheads bordering the triangular section are investigated and are 
shown to be important. 


INTRODUCTION 


Experimental investigations of swept box beams (references 1 and 2) 
have shown that the stresses and distortions in a swept wing can be 
appreciably different in character from those that would exist if the 
root were normal to the wing axis. The principal effect of sweepback on 
the stresses occurs under bending loads and consists in a concentration 
of bending stress and vertical shear in the rear spar near the fuselage. 
With regard to distortions, the effect of sweep is to produce some twist 
under loads that would produce only bending of an unswept wing and some 
bending under loads that would produce only twist of an unswept wing. 
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A relatively small amount of theoretical work has thus far been 
published on the analysis of swept wings of the shell type. Reference 3 
presents an energy method for determining influence coefficients of a 
swept box beam with bulkheads parallel to the flight path, and refer- 
ences 1 and 2 contain a semiempirical method for the deflection analysis 
of a swept box beam with bulkheads perpendicular to the spars. No pub- 
lications are known to be available on the stress analysis of a swept 
box beam with either type of bulkhead arrangement. 

The purpose of this paper is to present a method for the calcula- 
tion of both stresses and distortions of a swept box beam representing 
the main structural component of a swept wing with a carry-through sec- 
tion and with bulkheads perpendicular to the spars. The analysis is 
based on the four-flange idealized structure illustrated in figure 1. 

In a four -flange box beam of this type the carry-through and outer sec- 
tions can be analyzed by existing methods for unswept four-flange box 
beams if some means is provided for analyzing the triangular section 
which connects them. 

A method is presented for analyzing an idealized representation of 
the triangular section and for establishing continuity between the parts 
of the box beam; consideration is also given to the relationships between 
the idealized and actual structures and a comparison between the stresses 
and distortions calculated by this method and the experimental data of 
references 1 and 2. In the discussion, the effects of shear lag, which 
the method cannot give, are considered and an extension of the basic 
approach to permit their inclusion is indicated; also, the importance to 
the analysis of including the shear and bending flexibility of the bulk- 
heads bordering the triangular region is demonstrated. A complete 
numerical example is worked out in an appendix. 


SYMBOLS 

Principal Concepts 


A 

a,a n 

a ij 

^n^^ n ,D n 

b 
b 1 


area, square inches 
length of bay, inches 
coefficients of matrix 

arbitrary constants in solution of a differential equation 

width of outer section, inches 

width of carry-through section, inches 
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c 

E 

F 

f,J,P 


G 

h 

I 

J 

k i 

L 

l 

M 

N,N' 

P 

Q 

<1 

R 

T 

t 

U 

u 

V 

w 

X 


depth of box beam or bulkhead, inches 
modulus of elasticity, psi 

force group statically equivalent to a bending moment, pounds 

warping stiffness parameters 

modulus of rigidity, psi 

width of plate, inches 

moment of inertia, inches^ 

torsion constant, inches^ 

stiffness factors of outer and carry -through sections 

length of outer section, inches 

length of beam, inches 

bending moment, inch -pounds 

summations used in appendix B 

axial load in flange or stringer, pounds 

area moment, inches^ 

shear flow, pounds per inch 

ratio which has the value +1 for symmetrical loads and -1 for 
anti symmetrical loads 

torque, inch -pounds 

sheet thickness, inches 

strain energy, inch-pounds 

displacement in the x direction, inches 

vertical shearing force, pounds 

downward displacement or deflection, inches 

self -equilibrating, statically indeterminate force group, pounds 
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x,y,z rectangular coordinates, inches 

y distance from front spar to a specified center, inches 

a, p, 8, e,u) stiffness factors of beams 
7 shear strain 

£,ti nondimens ional parameters used in discussion of idealization 

9 angle of twist, radians 

A angle of sweep, degrees 

X effective width, inches 

H Poisson's ratio (assumed to be l/3) 

v effectiveness factor 

cr normal stress, psi 

t shear stress, psi 

t 

$ stress function 

♦,0 joint rotations, radians (see fig. 4) 

Subscripts 

Subscripts are used chiefly to designate the location of a dimen- 
sion, force, or stress, as follows: 

b cover of the box beam 

c web of the box beam or bulkhead 

F,f front spar or flange (see fig. 7) 

L left-hand end of beam 

Z along length of beam 

Z value at x = Z ] 

> (appendix B) 

o value at x = 0 
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p plate 

R,r rear spar or flange (see fig. 7) 

R right-hand end of beam 

x,y coordinate axis 

s stringer 

1,2,3* • -n specific locations shown in figures^ also, numbers to identify 

stiffness factors 

The singly exception to the foregoing convention is: 

e effective when applied to area, thickness, or moment of 

inertia 

Superscripts 

Superscripts are used to designate stresses and distortions 
produced by different types of action, as follows: 

B bending 

F F -force group 

R rigid-body displacements 

T torsion 

W warping 

X X-force group 

o flexure 

t shear 

ANALYSIS OF THE IDEALIZED STRUCTURE 
Basic Assumptions 

The type of idealized structure analyzed is shown in figure 1. It 
is a four-flange box beam, which is considered swept back in order to 
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avoid ambiguity in the designation of the front and rear spars. The 
sweptback parts are made up of two triangular sections and two rectan- 
gular outer sections which are symmetrical about and continuous with a 
rectangular carry-through section representing the part of the wing 
within the fuselage. The outer and carry-through sections contain bulk- 
heads which are placed perpendicular to the spars but the triangular 
section contains no bulkheads. All cross sections are symmetrical about 
a horizontal plane through the mid-depth of the spars. 

The box beam is supported, either rigidly or elastically, at the 
four corners of the carry-through section so that the reactions are 
simple vertical forces. It is loaded by a series of vertical forces as 
shown in figure 1. The resulting stresses are within the elastic range. 

The longitudinal direct stress in the idealized structure is assumed 
to be carried only by the concentrated areas at the corners of the cross 
section, and the side walls (spar webs) and covers are assumed to support 
shear stress only. The shear flow in the triangular cover sheets is 
assumed to be constant throughout the element (q^ in fig. 2). This 
assumption implies the existence of uniformly distributed normal forces 
on the hypotenuse of the triangular cover sheet but, since the adjacent 
carry-through section can carry normal force only at its corners, these 
distributed forces are lumped into two statically equivalent concen- 
trated forces (P 5 in fig. 2) at the ends of the hypotenuse. This 
assumption of uniform shear stress in the triangular sheet is approxi- 
mately justified by the experimental data of references 1 and 2 . 

The two bulkheads which border the triangular section are assumed 
to be beams with finite shear and bending stiffness in their own plane 
but with no resistance to distortion out of their plane. 

The relationship between the idealized structure described above 
and an actual structure is discussed subsequently. 


Method of Analysis 

The method of analysis is based on the assumption that the outer 
and carry-through sections can be analyzed by existing methods and that 
all that is then required is a means of establishing continuity between 
them through the triangular section. In order to accomplish this result, 
the structure is divided into a number of component parts as shown in 
figure 2. The forces assumed to exist on the cut sections are also 
shown. 

The first step in the analysis is to consider the vertical edges of 
the triangular section (joints 1 , 2 , and 3 ) as free bodies (fig. 3 ) and 
to write equilibrium equations for them. These equations include two 
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for moment equilibrium and one for vertical shear equilibrium at each 
joint or a total of nine equations. The number of equations is con- 
siderably less than the number of unknown forces , and the problem is 
therefore statically indeterminate. The number of equations, however, 
is exactly equal to the number of displacements required to specify the 
attitude and position of the three joints. These displacements include 
two rotations and a vertical translation of each joint, as shown in 
figure k. 

The next step in the analysis, therefore, is to establish force- 
displacement relationships for each component of the structure shown in 
figure 2. Through these relationships the forces appearing in the 
equilibrium equations can be replaced by the joint displacements and the 
loads applied to the structure. The nine equilibrium equations then 
contain as unknowns only the nine joint displacements and can be solved 
simultaneously for the displacements. Once the joint displacements are 
known, the force -displacement relationships can be used again to deter- 
mine the stresses and distortions of the entire structure. 

If one, or both, of the bulkheads (1-3 or 2-3) is .assumed rigid 
in its own plane, certain relationships among the joint displacements 
are immediately evident j thus the number of equilibrium equations needed 
is reduced and the analysis is simplified. 

In an analysis of this type many of the factors involved depend 
upon the nature of the applied load (symmetrical or anti symmetrical, 
bending or torsion) and it may therefore be advantageous to make a 
separate analysis for each type of load and then superimpose the results 
to obtain the desired solution. For convenience in the detailed develop- 
ment which follows, however, provisions for both bending and torsion 
are included simultaneously but with restrictions that they are either 
symmetrical or anti symmetrical about the carry-through section. 


If the three joints shown in figure 3 are considered as free bodies, 
a total of nine equilibrium equations can be written, two for moments 
and one for vertical shear at each joint, as follows: 


Joint-Equilibrium Equations 


Joint 1: 


Py - - Pj cos A = 0 


(1) 


p 9 + P 1 sin A = 0 

v l + C (<16 * <19 - 'll) = 0 


( 2 ) 


( 3 ) 
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Joint 2: 


Joint 3: 


p 2 " p 10 " 0 

w 

'tJ 

H 

H 

II 

O 

(5) 

V 2 + c( 4l - % - q 10 ) = 0 

( 6 ) 

Pq + P 4 cos A - P 3 sin A = 0 

(7) 

Pg - Pcj - Pjj. sin A - P 3 cos A = 0 

( 8 ) 

v 3 + c + iio " *8 + * 3 ) = 0 

(9) 


Inasmuch as the number of unknown forces appearing is greater than 
the number of equations, the problem is statically indeterminate. The 
principle of consistent displacements will therefore be used to obtain 
a solution. 


Force-Displacement Relations 

The attitude and position of the joints can be completely described 
by nine joint displacements, two rotations, and a vertical translation 
at each joint (fig. k) . Thus the number of unknown joint displacements 
is exactly equal to the number of equilibrium equations, so that a solu- 
tion is possible if sufficient force -displacement relations can be writ- 
ten to express all the unknown forces in terms of the nine joint 
displacements. 

All the internal forces (P's and q's) shown in figure 2 can be 
expressed in terms of the unknown joint displacements and the loads 
applied to the outer section, with the use of the force-displacement 
relationships for each component of the structure . All that remain axe 
the three vertical forces (V's) at the joints which are a combination of 
the vertical loads applied to the triangular section and the statically 
indeterminate support reactions. Since these forces are dependent upon 
the nature of each individual problem, they will be temporarily treated 
as known quantities) the modifications required for different types of 
supports are discussed in a subsequent section. 


The force -displacement relationships for each component can be 
written as indicated in the following sections. 
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Beams . - The two bulkheads 1-3 and 2-3 and that part of the front 
spar 1-2 bordering the triangular section can be analyzed as beams sub- 
jected to end shears and moments plus a running shear along the flanges. 
This running shear results from the shear flows in the covers adjacent 
to the flanges. The loading and distortion of a beam of this type is 
illustrated in figure 5» In appendix A , this type of beam is analyzed 
and the following general expressions are obtained for the end loads in 
terms of the end displacements, the running shear, and certain stiffness 
parameters which include both the shear and bending resistance of the 
beam: 



( 10 ) 


Wt — Wtj 

P L = -O0 L + P0 R - (a - (3 )~— r R - 5q z 


(ID 


P R = -00L + + (a - p) 


wl - wr 

— 1 — 


+ 6q z 


( 12 ) 


where 





€ 




and 


El 


bending stiffness of beam 


Get 


shear stiffness of beam 


§ 
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Specific force-displacement relationships for each of the beams can 
be obtained by application of these general equations with the following 
results: 

Front spar 1-2: 


q x = (a - 3) 12 


cos A + sin A + 0 2 
b tan A 


+ 2(a - p) 12 


vi - V2 

b 2 tan 2 A + fl2q 5 


(13) 

P 1 = - a 12(tL cos A + 02. sin AJ + 3 12 0 2 - (a - 3)12 23 tan A + 5 12 q 5 

^2 = "^12(^1 COS ^ + ^1 s ^- n + a 12^2 + * P)l2 h tan a" ~ 6 12 q 5 (15) 

Bulkhead 1-3: 

0-i - \|ro cos A + 0o sin A v-i - Vo 

*9 - <* - P)i3 + 2 <“ - P)i 3 ~ 

e l3 ( q 7 + q ll) 

P 9 = - 0,1301 - 313 (*3 cos A - 03 sin a) - (a - 3)13 
5 13 ( q T + q ll) 

/ , \ Vi - Vo 

Pq = -^ 2 . 3^1 * <x 13 (^3 cos A ~ ^3 sin A j + ( a - P)j- " — - — - + 


'13 b sec A 

s i 3 ( 9 7 + ’ll) 

Bulkhead 2-3: 

l|f 2 + 1^3 Wo — w 2 

q 10 = “ (a - P) 2 3 ~ 2(a - P) 2 3 + e 23( q 4 - q 5) 

P 4 = - a 23*23 + *23*2 - (a - p) 2 3 “ 5 2 3 ( q 4 - %) 

P 11 = '^23*3 + a 23^2 + ( a " ^23 ""^T> + 5 23 ( q 4 ■ *%) 


(16) 


(17) 


( 18 ) 

(19) 

( 20 ) 

( 21 ) 
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In these equations the subscripts 12, 13, and 23 applied to the 
stiffness parameters a, 3, &, and e refer to the front spar 1-2 

and the two bulkheads 1-3 and 2-3, respectively. Several unknown forces 
namely, qj^, q^, q^, and q-Q, which do not occur in the joint equi- 

librium equations, appear on the right-hand sides of equations ( 13 ) to 
(21) j these forces, too, can be expressed in terms of the nine joint dis 
placements when the other structural components are considered. 

Triangular cover sheet . - The triangular cover sheet is assumed to 
carry a uniform shear flow q^ along its mutually perpendicular edges 

(1-2 and 2-3). In order that this element be in equilibrium, shear and 
normal stresses are required along the hypotenuse and the corresponding 
forces are shown (fig. 2) as a uniform shear flow q-Q acting along 


that edge and a pair of concentrated forces P 5 
The equilibrium equations are: 

acting at the joints. 

qil = -qcj cos 2A 

( 22 ) 

P^ = -q^b sin A 

( 23 ) 


Force -displacement relationships are obtained by assuming that the 
maximum shear strain in the sheet is equal to the amount by which the 
right angle 1-2-3 is changed. In terms of the joint rotations, this 
shear strain is: 

clh sln A - 01 cos A - * 2 0 3 - 0 2 \ /rtllX 

7 l-2-3 2 \ b tan A b / ' ' 

Then, 

<35 = cos A - cot A - cos A cot A - ^ + (25) 

and the relations for q^ and P follow Im mediately from equa- 
tions ( 22 ) and ( 23 ). 

Outer section .- That part of the structure outboard of bulkhead 2-3 
acts as an unswept cantilever box beam supported on a flexible root and, 
as such, can be analyzed by existing methods of analysis. The stresses 
and distortions at any point can be expressed in terms of the applied 
loads, the distortions of the root, and certain elastic stiffness 
factors. Then, the force -displacement relationships required to define 
the internal forces at the root are: 
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( 30 ) 


( 29 ) 


( 27 ) 


(28) 


(26) 


In these equations V, M, and T represent, respectively, the applied 
vertical shear, bending moment, and torque (about some reference axis) 
at the root of the outer section and the k's represent elastic stiff- 
nesses of the outer section. The stiffness factors kq, k2, k4 and 
the like are functions of the distribution of the applied loads and the 
dimensions and material of the outer section, whereas k 3 and ky 
depend only on the latter . The quantity 02 - 03 is a measure of the 
warping of the root cross section and is the only root distortion 
appearing in the equations, since the others are rigid-body movements 
which do not affect the stress distribution. Thus, effectively, the 
root bulkhead is assumed rigid in its own plane as far as the outer- 
section analysis is concerned. 

Any method of analysis can be used to determine the stiffness 
factors provided that cross-sectional warping and its restraint are 
taken into account. This provision requires a more refined approach 
than is made in elementary bending theory. The stiffness factors are 
the same for symmetrical and antisymmetrical loadings but, since bending 
and torsion produce different types of effects, they have been separated 
in the equations. In order to evaluate the torque T, the loads must be 
referred to a reference axis. The most desirable axis is one which makes 
the stresses at the root due to the bending moment M equivalent to 
those given by elementary theory, although it is not generally possible 
to achieve this relationship at all stations. The so-called "shear 
center does not locate such an axis. The choice of a reference axis 
will be treated at greater length in the section on idealization. 

Carry-through section .- The carry-through section, like the outer 
section, is a box beam that can be analyzed by existing methods. In 
this case, however, the stress distribution is expressed in terms of 
only the end distortions since internal end forces are the only loads 
applied. The force -displacement relationships are then: 

Pg = + ijf^ sin A + 0^ cos A j + k -| p ^ sin A - 0^ cos A^ + 

kl 3 ^wq + W 3 ^ + kqt(wq - W 3 ) 


( 31 ). 
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P'7 = 


% 


q 7 


k l5 (^1 + ^3 sin A + ^3 cos A ) + k i6(^i " ^3 sin A - 03 c° s A ) + 
k l7 ( W 1 + ^3) + k l8( w l " w 3) (32) 

^19(^1 + ^3 sin A + 03 cos A ) + ^0(^1 “ ^3 sin A " 03 cos A ) + 
k 21 ^ w i + w 3^ + k 22( w l " w 3) (33) 

^23 (^1 + ^3 s ^ n A + 03 cos A )+ ^4(^1 “ ^3 s ^ n A - 03 cos A ) + 

^ 25 ( W 1 + w 3) + ^6^1 " * 3 ) (34) 


^ sin A + 0^ cos a) + k 2 g(ii “ ^3 sin A - 03 cos A ) + 

* 3 ) + k 3 o( w l “ v 3 ) (35) 

In these equations, the k f s represent elastic stiffness factors, which 
depend upbn whether the loading is symmetrical or antisymmetrical, as 
well as upon the dimensions and material of the carry- through section. 
They may he determined by any method of analysis as long as cross- 
sectional warping and the shear and bending stiffness of the spars are 
considered. The displacements which appear have been so grouped that 
they have a particular physical meaning. Thus, the quantities associated 
with the first, second, third, and fourth terms in parentheses represent 
a bending type of rotation, a warping, a translation, and a torsion type 
of rotation of the end cross section, respectively. 



Solving for the Joint Displacements 

The force -displacement relationships (equations (13) to (35)) are 
sufficient to express all the internal forces in the equilibrium equa- 
tions in terms of the nine basic joint displacements and the applied 
loads. Upon substitution, the nine equations involve only nine unknown 

W, Wp w^ 

displacements; namely, \|r 2 > + 3 , 0 1 , 03 > and ^ ' 

and they can then be solved numerically for these unknowns. 

The equations obtained by direct substitution have coefficients 
containing many terms which are tedious to evaluate; however, a number 
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of combinations can be made which substantially simplify the final 
equations. Equations (l) to (9) are combined as follows to obtain nine 
simpler equations; 


[(l) - (8)J sec A - (4) + ( 5 ) tan A 

(36) 

[(2) + (T)J esc A + (4) + ( 5 ) cot A 

(37) 

§[<3) + ( 6 ) + ( 9 )\ 

(38) 

(4) 

(39) 

(5hp- 

°23 

( 40 ) 

— ( 6 ) tan A 

(ia) 

(7) sec A + (5) 

( 42 ) 

( 8 ) sec A - ( 5 ) tan A 

(43) 

^(9) sec A 

(44) 


The resulting system of equations is written in matrix form as 
follows: 


a ll 

a 12 

a l3 

a l4 

a 15 

a l6 

a 17 

a l8 

a 19 


*1 


a 10 

a 21 

a 22 

a 23 

£*24 

a 25 

a 26 

a 27 

^8 

a 29 


*2 


a 20 

a 31 

a 32 

a 33 

a 34 

a 35 

a 36 

a 37 

a 38 

a 39 


^3 


a 30 

a 4l 

a 42 

a 43 

a 44 

a 45 

a 46 

a 47 

a 48 

a 49 


01 


a 4o 

a 51 

a 52 

a 53 

a 54 

a 55 

a 56 

a 57 

a 58 

a 59 


02 

= 

a 50 

a 6l 

a 62 

a 63 

a 64 

a 65 

a 66 

a 67 

a 68 

a 69 


03 


a 6o 

a 71 

a 72 

a 73 

a .74 

a 75 

a 76 

a 77 

a 78 

a 79 


wi 

b 


a 70 

a 8l 

a 82 

a 83 

a 84 

a 85 

a 86 

a 87 

a 88 

a 89 


W2 

b 


a 8o 

a 9i 

a 92 

a 93 

&94 

a 95 

&96 

a 97 

a 98 

a 99 


w^ 

b 


a 90 

u_ 


The coefficients aij are given by the expressions in table I. 
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The terms in the matrix of coefficients involve the elastic stiff- 
ness factors and dimensions of the structurej the constant terms contain 
elastic stiffness factors, the loads applied to the outer section, and 
the vertical forces at the joints. Each of these joint forces contains 
a component of the load applied to the triangular section and, in addi- 
tion, Vq and V 3 contain the support reactions which may he statically 
indeterminate. The loads applied to the triangular section are so 
divided among the three vertical edges that the resulting forces form a 
statically equivalent system. The reactions depend upon the nature of 
the supports and are introduced into the analysis as indicated in the 
following sections. 

Rigid supports .- In the case of rigid supports wq and wq are 
zero and there are thus seven unknown joint displacements which require 
only seven equations for their determination. In the matrix, equa- 
tion (45), columns 7 and 9, which are the coefficients of wq and W3, 
respectively, can he immediately eliminated. The required seven equa- 
tions are then obtained hy the elimination of two rows, the logical ones 
being rows 3 and 9 since they are derived from equations containing the 
unknown support reactions. After the joint displacements have been cal- 
culated, Vq and V3 can be determined by substitution into equa- 
tions ( 3 ) and ( 9 ). 

Elastic supports . - For the case of elastic supports, all nine equa- 
tions ar 2 required but must be modified to include force -displacement 
relationships for the supports. The joint forces can be expressed as 
follows: 


Vi - k 31 V! + Vi’ 

V 3 . k 32 w 2 + Vj' 

where 

k support stiffness factor 

V 1 component of loads applied to triangular section 


Calculating Stresses and Distortions Throughout 
the Idealized Structure 

The stress and distortion distributions for the complete box beam 
have been defined in terms of the applied loads and the nine joint dis- 
placements. Once these joint displacements have been determined by 
solving equation (45), the procedures outlined previously can be reversed 
and all of the forces at the joints can be calculated. 
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Determination of the detailed distributions is slightly more com- 
plicated. For the front spar and the two bulkheads bordering the tri- 
angular section, the equations of appendix A can be used. For the outer 
and carry-through sections, the complete stress distribution can be 
determined from the analysis that was used to obtain the stiffness 
factors kp to k30- The effects of rigid body motions, which do not 
affect the stresses, must be included in the calculation of distortions. 

The relationships between the computed stresses in the idealized 
structure and the actual structure are discussed in the section on 
idealization, which follows. 


IDEALIZATION OF AN ACTUAL STRUCTURE 
Outer and Carry -Through Sections 


The outer and carry-through sections are unswept box beams which 
can be analyzed by existing methods of analysis. Since such methods are 
by no means standard, however, a definite procedure is presented in 
order that the idealization of the complete structure may follow a 
consistent pattern. 

The basic assumption that the idealized outer and carry-through 
sections are conventional four-flange boxes implies that the normal 
stress in the walls of the actual box beam varies linearly between 
adjacent corners. A generalized stress distribution of this type can 
be represented by a linear combination of the two stress distributions 
shown in figure 6(a), one of which equilibrates the applied load and is 
uniform across the cover, while the other is self -equilibrating and 
varies linearly across the cover. The uniform distribution is designated 
bending stress because it is obtained from elementary beam theory which 
assumes that plane cross sections remain plane after loading. Similarly, 
the linearly varying stress is designated warping stress because it is 
associated with the warping of the cross section out of its plane. 

The normal stresses on the actual cross section are represented by 
four concentrated forces at the corners of the idealized cross section. 
The total force at each corner consists of two components, one from the 
F-force group corresponding to the bending-stress distribution and one 
from the X-force group corresponding to the warping-stress distribution, 
as shown in figure 6(b). The equivalence between the force group and 
the corresponding distributed stress is determined on the basis of over- 
all statics of the cross section and the moment applied to each cover. 

The effective flange areas of the idealized structure axe then chosen so 
that the flange stress in the idealized structure is equal to the corner 
stress in the actual structure. 
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A cross section of the type shown in figure 7(a) can then be 
idealized as follows: 

(l) Obtain the equivalent cover (fig. 7(b)) by adding, to the 
actual cover and corner flanges, areas representing the moment -carrying 

capacity of the webs; that is, ^ct r and ^ctf; the area of the 

equivalent cover is therefore: 


Ab = Af + ictf + A r + ict r + btb + 1 .... (As). 


n 


6 T T 6 V 

(2) Locate the centroid of the equivalent cover: 


(46) 


? = + Ar + Z ct r + ZT_ 


n 


As b 


n 


(Vf) 


(3) Calculate the moment of inertia of the equivalent cover about 
a vertical axis through its centroid: 

^ - ( Af + 5 ctf )(lf * ( Ar + - if * # tb * btb (§ - If + 

E A s n (v - if <W) 


(4) The effective area of each front flange is then: 
(a) For bending stresses: 


- (1 - IK - ( b - 1)1 <*> 

where I is the moment of inertia of the entire cross section about 
the horizontal axis of symmetry. 

(b) For warping stresses: 


W _ Ib/ft 2 
y/b 


( 50 ) 
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(5) The effective area of the rear flange in each case is: 


4 



\ 



(51) 


If many equally spaced stringers are used, satisfactory results 
can be obtained by treating them as an equivalent sheet and thus 
eliminating the evaluation of lengthy summations. 


It is important to note that a different effective area is 
associated with each type of stress distribution, as should be expected, 
since each is associated with a different type of physical action; 
therefore, if accurate results are to be obtained, the two types of 
stress distributions must be completely separable in the analysis, that 
is, they do not appear simultaneously in the evaluation of any one 
stiffness factor in equations (2 6 ) to (35)- This separation is not 
generally possible; however, one way to accomplish complete separation 
in the outer section will be described. Similar considerations apply 
to the carry-through section. 


The outer section is an unswept cantilever box beam on a flexible 
root and the forces on any cross-section as given in equations ( 26 ) 
to ( 30 ) can be expressed as the sum of: (l) forces that exist in the 

loaded cantilever on a rigid root and (2) forces that exist in an 
unloaded cantilever having the root warped an amount (02 “ 03)- Since 

root warping produces only warping stresses, the effective areas for 
warping stresses (equations ( 50 ) and ( 51 )) are used for the determina- 
tion of the stiffness factors k 3 and ky. The choice of effective 
areas for the analysis of the loaded cantilever is more difficult because 
the application of vertical loads will, in general, produce both bending 
and warping stresses; however, since torque loads produce only warping 
stresses, it may be possible to locate some axis along which applied 
vertical loads will produce only bending stresses at every cross section; 
then, the loading can be divided into vertical forces applied along this 
axis and torques about this axis. The stress types are thus separated 
and the effective areas for bending stress (equations (*4-9) and ( 51 )) can 
be used to calculate the stiffness factors associated with the vertical 
loads (kj_, k^, k^, kg and k 10 j and the effective areas for warping 
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stresses (equations ( 50 ) and ( 51 )) can be used to calculate the stiff- 
ness factors ^kg, kg, and associated with the torque. 

An axis of the type described obviously exists for a doubly 
symmetrical cross section; such an axis is also known to exist for a 
four -flange box beam of constant cross section which is symmetrical 
about a horizontal plane (fig. 7(c)). The location of this axis, at the 
center which might be called the zero warping center of the cross 
section, is given by: 

£ = < 52 > 

where 

1 _ = i/l -L' \ 

t c 2 \tf + t r j 



lb tc 


The zero warping center should not be confused with the shear center. 
Vertical loads applied along an axis through the shear center will 
deflect the box beam without twisting it, if the cross sections are free 
to warp, a condition that is not satisfied at the rigid root of a 
cantilever. Vertical loads applied along an axis through the zero 
warping center will result in a combination of deflection and twist, but 
the cross sections of the box will not warp. If the cross-section is 
doubly symmetrical, the zero warping center and shear center coincide at 
the geometrical center of the cross section. 

The preceding discussion has been exclusively concerned with the 
problem of converting the actual structure into an idealized structure 
that can be easily analyzed. After the analysis has been completed and 
the magnitudes of the corner forces determined, the problem of con- 
verting corner forces into stress distributions arises. This conversion 
of forces is accomplished by determining the stresses corresponding to 
each type of force group (F or X) and then summing the stresses to get 
the total stress. The relationship between force groups and stress 
distributions are shown in figure 6. The type of force group is deter- 
mined from physical considerations; for example, equation ( 26 ) shows 
that P 3 is composed of three forces of which k]_M is of the bending 
type whereas k^T and k^ ($2 - 0^) are the w &rping type. 
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With regard to shear stresses, the analysis of the idealized 
structure gives the average value of the shear flows in the walls of 
the actual structure. A more detailed shear-stress distribution can be 
derived from the distributed normal stresses; however, this additional 
refinement is probably unwarranted in view of the many approximations in 
the basic solution. 


Triangular Section 

The idealized triangular section consists of three parts: the two 

cover sheets which are in a state of uniform shear, and the front spar, 
which is assumed to be a beam. 

The thickness assigned to the idealized cover sheet should properly 
represent the shear stiffness of the actual cover. For unbuckled sheet 
alone, this thickness is that of the actual sheet. This value should be 
decreased if the sheet has buckled or increased if there are closed- 
section stringers which contribute to its shear resistance. Similar 
considerations apply to the determination of the shear resistance of 
any other element of the structure. 

In the calculation of the moment of inertia of the front spar, a 
contribution from the sheet and stringers in the cover of the triangular 
section must be included to account for their ability to carry direct 
stress. For simplicity, the moment of inertia of the idealized front 
spar is assumed constant in the spanwise direction and thus may be deter 
mined by treating the triangular bay as a rectangular bay of constant 
cross section equal in width to the triangular bay where it joins the 
outer section; an effective area may then be assigned to the idealized 
front spar by the method recommended for the outer section. This ideal- 
ization will result in two moments of inertia, one for bending and one 
for warping stress. Again, separation of these two stress systems in 
the analysis is desirable but in this case it is only partly possible. 

If the load on the outer section is torsion only, only warping stress 
will exist in the triangular section and the effective warping area 
should be used. If the loads on the outer section are of the bending 
type only, both kinds of stresses will exist in the triangular section 
and direct separation is impossible; however, the warping stress is 
usually small compared with the bending stress and in such cases 
satisfactory results can be obtained by using the effective bending area 


Bulkheads 

When values of flexural stiffness are assigned to the two bulkheads 
bordering the triangular section, consideration must be given to the 
fact that bending of these beams is accompanied by extension or 
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compression of the adjacent cover sheets. The analysis in appendix B 
gives the distortions of a combined bulkhead and cover sheet and the 
moment of inertia of a simple bulkhead which has the same distortion. 

The results can be summarized as follows: 

I e = 1(1 + V) (53) 


where 

I 

^e 

v 

c 

7 

^e 

E 


moment of inertia of bulkhead alone , inches 4 
effective moment of inertia of combination, inches^ - 


an effectiveness factor plotted in figure 8 as a function of 

,'c 2 Zt e ' 

the nondimens ional parameter 


depth of bulkhead, inches 
length of bulkhead, inches 

equivalent thickness of the cover sheet, inches 


E sheet 

E bulkheady 


modulus of elasticity, psi 


COMPARISON BETWEEN THEORY AND EXPERIMENT 


The accuracy of the method is demonstrated by comparing calculated 
stresses and distortions with the test data of references 1 and 2 . The 
test specimen used is illustrated in figures 9 and 10 and the details 
of the calculations are given in the numerical example of appendix C. 

The comparisons are presented graphically for each of the four test 
conditions, symmetrical and antisymmetrical tip bending and torsion 
loads as follows: 

(a) Distortions of the outer section in figure 11 

(b) Spanwise distribution of spar shear stresses shown in the left- 

hand parts of figure 12 

(c) Spanwise distribution of flange normal stresses shown in the 

right-hand parts of figure 12 

(d) Chordwise distributions of normal stress at three spanwise 

stations in figure 13 
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In each case the sign conventions employed are those of references 1 
and 2 which are occasionally in conflict with those employed elsewhere 
in this paper. 

The test data in these figures are presented in the usual manner 
and several calculated curves are given to illustrate different theo- 
retical approaches, as follows: 

(a) All figures contain a heavy solid line which represents the 
results of the numerical example of appendix C. The area under this 
curve is vertically hatched in the stress plots (figs. 12 and 13 ). 

(h) A dash-dot line appears in some figures to show the effect of 
superimposing shear-lag effects on the results of the numerical example. 
The determination of these effects is described in the discussion. 

(c) The stress plots also contain dash lines which give results 


In general, the results of the numerical example (solid lines) are 
in fair agreement with the test dataj however, much better agreement is 
achieved when shear-lag effects are added where applicable. Elementary 
theory gives the least satisfactory results since it does not include 
the effects of either sweep or shear lag. The discrepancies between the 
calculations of the numerical example and the experimental data are 
primarily the result of analyzing an overly simplified idealization of 
the actual structure. The assumed idealized structure is incapable of 
distorting in all of the shapes assumed by the actual box beam; therefore, 
the analysis cannot give completely accurate results. The most signifi- 
cant effect of oversimplification is the neglect of shear-lag stresses. 
Shear lag appears whenever the webs carry vertical shear stress and it is 
characterized by an increase in normal stress in the vicinity of the 
flanges with a corresponding decrease in the rest of the cover (see 
fig. Ik). There is also an associated change in the shear stresses in 
the cover. 

The effect of shear lag on the stress distribution is most evident 
in the chordwise plots of figure 13 . It appears in the outer section 
for the bending loads and in the carry-through section for the anti- 
symmetrical loads. In each case the effect is carried over into the 
triangular section because of continuity. 

Shear-lag stresses effectively reduce the stiffness of the structure 
and thus increase its deflection. The cantilever beam deflection of the 
outer section is increased by its own shear lag, but for this specimen 
the increase was small enough to be neglected. More important is the 
reduced stiffness of the carry -through section under anti symmetrical 


obtained from elementary theory 
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loads which results in an increased rotation of the ends of the carry- 
through section. This effect causes only a rigid-body rotation of the 
triangular and outer sections since the carry-through section is doubly 
symmetrical. The effect of shear lag in the carry -through section was 
estimated by reducing the moment of inertia of the carry-through section 
(l e = 0.2521) by the procedure described subsequently in the ’’Discussion" 
with the results shown in figure ll(b) and 11(d). For the antisymmetri- 
cal bending load the effect was underestimated whereas it was over- 
estimated for the anti symmetrical torsion load. 

Some other effects of oversimplification are associated with the 
idealization of the triangular section. The actual structure had a 
short bulkhead (bulkhead 7> fig. 9) in the triangular section which was 
neglected in the analysis. Its presence introduces additional restraints 
which change the shear stress in the front spar (fig. 12) and the chord- 
wise normal-stress distributions in the carry-through section (fig. 13 ). 
Also, a number of approximations were used with regard to the effective 
moment of inertia of the front spar which introduce uncertainties in 
the analysis. 

Another factor that affects the agreement is the flexibility of 
the supporting jig used in the tests. This jig was assumed rigid in 
the analysis but deflected during the tests. Under symmetrical loads 
(reference 1) these deflections amounted to a rigid-body rotation of 
the complete structure and it was a simple matter to correct the meas- 
ured deflections; however, under antisymmetrical loads, a small amount 
of twist remained in the carry-through section after the rigid-body 
corrections had been made (reference 2). The method of analysis devel- 
oped in this paper was used to calculate the effect of the measured 
carry-through section twist on the theoretical stresses and deflections. 
For the antisymmetrical bending load the principal changes were in the 
shear stresses in the carry-through section, as might be expected, 
because of the torque required to twist this section. The warping 
stresses also changed throughout the box beam and the deflection of the 
outer section increased. For the antisymmetrical torsion load the twist 
was small enough to be negligible. In general, these changes improved 
the agreement between theory and experiment but were not of sufficient 
magnitude to warrant their addition to the calculated results. 


DISCUSSION 

Determination of Shear -Lag Effects 


The method presented in this paper is for the analysis of a four- 
flange box beam which experiences only first-order warping. Thus, the 
stress distribution in the actual structure varies linearly between 
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corners and consists of the bending and warping stress components shown 
in figure 6. An actual structure experiences other types of cross- 
sectional warping, a common one being the second-order or shear-lag 
type which introduces a departure from linearity in the cover stresses 
by superimposing a self -equilibrating stress distribution of the type 
shown in figure 14(a). The experimental data of references 1 and 2 
show that shear lag is important in the specimen tested because it 
causes an increase in the flange normal stresses which is, in some cases, 
larger than the changes caused by first-order warping and increases the 
flexibility of the structure which increases the deflections. Some 
means for calculating these effects is therefore necessary. 

Any method used to calculate shear-lag effects requires the 
analysis of a more complicated idealized structure than the four -flange 
box beam. Conforming to the previous assumptions regarding the stress- 
carrying ability of the idealized structure, the minimum addition is a 
single central stringer in each cover as shown in figure 14(b); other 
additional cover stringers permit calculation of the effects of third 
and higher order warping. These more complicated structures can be 
introduced into the analysis in either of two ways: (1) The basic 

method can be extended to the direct analysis of the more complicated 
structure and thus automatically include shear- lag effects. (2) The 
results of the simpler analysis can be corrected by a process which 
combines experimental data and individual shear- lag analyses of the 
outer and carry-through section. Each of these approaches is briefly 
described. 

The direct extension of the basic method involves the analysis of 
an idealized structure of the type shown in figure 15 . The analysis 
follows the procedure previously described but two new features are 
introduced. First, more joints are involved, for which additional 
equilibrium equations are required and thus a larger system of simul- 
taneous equations must be solved. Second, force -displacement relation- 
ships for the outer, triangular, and carry-through sections must be 
modified to account for the new types of forces and distortions of the 
idealized structure. 

The correction process can take on a variety of forms, two of which 
were used to calculate the shear -lag corrections applied to the results 
of the numerical example, appendix C, to obtain the dash-dot lines in 
figures 11, 12, and 13. The shear-lag corrections for the outer section 
were determined by using the single-substitute stringer method (refer- 
ence 4) to calculate the shear-lag stresses in the outer section. The 
outer section was analyzed as an ordinary, unswept box beam on a rigid 
root and a constant empirical multiplying factor was used to obtain 
good over-all agreement with the experimental spanwise and chordwise 
normal-stress distributions. The multiplying factor accounts for the 
root restraint provided the outer section by the triangular section; 
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this factor was found to he 0.77 for the symmetrical bending load and 
1.28 for the antisymmetrical bending load. 

A somewhat different procedure was used to determine the shear -lag 
corrections for the carry-through section because it was found that the 
shear-lag stresses, caused by the antisymmetrical bending load, could be 
approximated by a shear-lag analysis which assumed that the end cross 
sections were restrained from warping. This approach, however, was less 
satisfactory for the antisymmetrical torque load. The method used 
analyzed the cover as an equivalent sheet which carried both shear and 
normal stress (reference 5 ) and it gave better chordwise and spanwise 
stress distributions than a similar analysis which used the single sub- 
stitute stringer method. The reduced moment of inertia used to deter- 
mine the effect of shear lag on the distortions of the carry-through 
section was also obtained from this analysis. 

From considerations of accuracy, the preferred method for the 
determination of shear -lag effects is the direct extension to a more 
complicated idealized structure} however, such an analysis requires a 
large amount of work. The ease with which the correction process can 
be used is a definite advantage, but it can be applied, with assurance 
of accuracy, only to structures closely resembling the test specimen 
from which the empirical factors were obtained. Even then, the correc- 
tion process is only fairly accurate because it cannot adequately 
account for the interaction between the various parts of the structure. 


Effects of Bulkhead Flexibility 

In the analysis of unswept box beams the internal bulkheads are 
often assumed to be rigid in their own plane. This assumption yields 
satisfactory results except when a discontinuity of structure or loading, 
such as a cut-out, introduces large loads into a bulkhead. A study of 
the test results presented in references 1 and 2 shows that bulkheads 6 
and 8 (fig. 9) of the test specimen were subjected to substantial shear 
and bending loads; thus, their distortions may have an important effect 
on the structure. The shear and bending flexibility of the bulkheads 
bordering the triangular section is included as a basic feature of the 
method of analysis presented in this paper although the development 
could have been considerably simplified by assuming them rigid. 

The effect of bulkhead flexibility on the stress and distortions of 
the swept box beam of references 1 and 2 was investigated by solving a 
series of numerical examples similar to that of appendix C. These 
examples were for symmetrical tip bending and torque loads for the four 
cases of bulkhead flexibility listed in the following table (all bulk- 
heads other than bulkheads 1-3 and 2-3 (fig. 2) were assumed rigid): 
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Case 

Bulkhead 1-3 

Bulkhead 2-3 

I 

Flexible 

Flexible 

II 

Flexible 

Rigid 

III 

Rigid 

Flexible 

IV 

Rigid 

Rigid 


Selected stresses and displacements calculated for each case are 
compared in figure 16. Despite the fact that the bulkheads were con- 
structed of i-inch steel plate, appreciable errors occurred when they 

were assumed rigid in their own plane. The effects of bulkhead flexi- 
bility on the stresses were more pronounced for torsion than for bending 
loads} however, for each type of load the solution for case IV over- 
estimated the warping stress in the outer section. The different root 
distortions calculated for each of the four cases lead to slightly dif- 
ferent deflections and somewhat greater differences between the rotations 
of the outer section. 

In general, the results indicate that the flexibility pf the bulk- 
heads bordering the triangular section has an important effect upon the 
stresses and distortions of a swept box beam and should be considered 
in the analysis if accurate results are to be obtained. 

In addition to the studies of bulkhead flexibility, some investiga- 
tions were made of the effect of the number of bulkheads in the tri- 
angular section. For example, a numerical analysis which used the 
assumption of closely spaced rigid bulkheads, often used in shell 
analysis, gave very erroneous results for the stress distribution in 
the triangular section. In all the cases considered, the number of 
bulkheads in the triangular section had only a small effect on the 
stresses outside of the triangular section. The experimental data in 
figures 12 and 13 illustrate the effect of an extra bulkhead on the 
stresses within the triangular section. 


Extension of the Method 

The method of analysis in the form presented, is not expected to 
be generally applicable to the precise analysis of all swept wings 
because of differences in structural arrangement and the degree of 
idealization assumed} however, the basic approach can be used in many 
situations . 

Extension of the method to other types of swept wings is straight 
forward if the bulkheads are perpendicular to the spars. The modifica- 
tions required to obtain more accurate stress distributions have been 
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indicated in the section on the determination of shear-lag effects. 
Similar procedures are required for the analysis of multispar box beams, 
that is, additional joints are created where the spars cross the 
bulkheads. 

Application of the basic approach to swept box beams in which the 
bulkheads are placed parallel to the flight direction introduces some 
new problems. In the first place, a swept wing of this type does not 
have an outer section which can be analyzed by existing methods. Exten- 
sion of the method of joints, which is a feature of the basic approach, 
to the entire structure substantially increases the complexity of the 
solution and some other approach may be more desirable. A second prob- 
lem is the establishment of force -displacement relationships which cor- 
rectly predict the physical behavior of the parallelogram- shaped cover 
sheets . 


CONCLUDING REMARKS 


A method has been described for the stress and distortion analysis 
of a swept box beam with a carry-through section and with bulkheads 
perpendicular to the spars. The method is based on a simple four-flange 
box type of idealized structure and permits an estimation of the first- 
order warping stresses that result from sweep but does not permit the 
evaluation of higher -order stresses such as shear lag. Agreement with 
experiment is therefore only fair; however, extension of the basic 
approach to permit more refined analyses, which include shear-lag 
effects and other structural arrangements such as multiple spars, has 
been indicated. 

The method assumes that the outer and carry-through sections are 
unswept box beams and thus can be analyzed by existing methods. Con- 
tinuity is established between them through the analysis of the tri- 
angular section. The analysis of the triangular section isolates the 
structural joints as free bodies and gives an equilibrium equa.tion for 
each degree of joint freedom. The joint forces are expressed in terms 
of joint displacements and a set of simultaneous linear equations, which 
completely defines the joint displacements, is thus obtained. 

The method takes into account the flexibility, both as regards 
shear and bending, of the bulkheads in and around the triangular section. 
The results of a numerical study have been presented to show that 
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appreciable errors can appear in the calculated stresses and deflections 
if the usual assumption of rigid bulkheads is used in this region. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., August 17 , 1950 
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APPENDIX A 


FORCE -DISPLACEMENT RELATIONSHIPS FOR BEAMS 


The beam analyzed is assumed to be of constant idealized cross 
section with a web which carries only shear stress and subjected to the 
loading shown in figure 5(a)* This loading consists of flange axial 
loads Pl and Pr at the ends of the beam, a constant running shear 
flow q^ applied to the flanges and a constant shear flow q c in the 
web. The distorted shape of the beam under load can be described by the 
end displacements w^ and Wp and the end rotations and 0p 

(fig. 5(b)). 

Consideration of a differential element of the beam (fig. 5(c)) 
yields the following equilibrium equation 


dP 


- q 7 - q c = 0 


dx n 

and the following relations between loads and distortion 


(Al) 


*c - <*(* * I 

p = El d£ 
c dx 


(A2) 

(A3) 


Substituting equations (A2) and (A3) in (Al) yields 

dw EI_ d 2 0 q 2 

dx Get ^2 + ^ + Gt 

and, since q c is constant, from equation (A2), 

d0 _ d 2 w 


dx 


dx c 


(A4) 


(A5) 


Equations (A4) and (A5) have solutions which can be expressed as 
follows: 


0 = Cp + CpX + CoX 2 


w = Cj_x + |c 2 x 2 + c 3^ x ^ 


2EIx \ 
Get 


J + c k + 


*1 

Gt 3 


(a6) 

(A7) 
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where Cp, Cg, C3, C4 are constants which can be determined from the 
boundary conditions. In this case, the boundary conditions are 
prescribed as follows: 


w x=o = W L 
w x=Z = W R 
^x=o - 
^x=Z = 0R 




(A8) 


These boundary conditions require that the constants C n have the 
following values: 


Cl - 0 L 

C2 = y(0R - 0 l) - C 3Z 


C 3 " -2Z 


C 4 ~ W L 


2 ( u r - ”l) - + 0; 

[)- §^z 

Z 2 2EI 

6 + Get 



(A9) 


y 


An expression for the load P can be obtained from the substitution 
of equations (A 6) and (A9) in equation (A3), as follows: 


P = 


a (sx 2 4EI 


*4 


2x + 4EI 
1 GctZ 2 , 


+ 2 


'2x 


' W R 


- w L q 
+ 77 


l 


clfl + 4 ei 
EI \3 GctZ 2 


(A10) 


Similarly, from equations (A2), (A6), (A7), and (A9), 

|(0R + 0l) + ^(wl - «s) - 


QLr = 


Z 2 

cZ/1 4-EI 

EI \3 + Gctz * 


(All) 


NACA TN 2232 


31 


The applied loads can now he expressed in terms of the distortions 
if the foregoing equations are evaluated at the boundaries. These rela- 
tions for the loads can be conveniently written as follows: 


P = - aJf L + P0r - (a - . 5q 


W L - W R 


Pr - “P0L + a 0R + ( a “ 3) i + 


(A12a) 

(A12b) 


where 


- (a - * 2(a - (D^jgS 


^Z 


(A12c) 


a 


1/4 4ep \ 

®\3 GctZ 2 / 


3 


1 f_2 + 4EI \ 

Gctl 2 / 



(JO 


cz/l 4EI \ 
EI \3 + GctZ 2 / 


It is often desirable to express the distortions in terms of the 
loads. These expressions can be obtained from a few simple operations 
on equations (A12). Addition of equations (A12a) and (A12b) gives the 
following relationship between end rotations: 


0R = 0L 

= 0L 


P L + P R 

4 - — — 

a 4 - 0 


iifr + %) 


(A13) 
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Substitution of (A13) into (A12a) or (A12b) yields the following 
alternate expressions for the difference in end displacements 


vr - wl = 0L* + + 

W R ' W L = 0R Z + 75T + L 




PL - PR . cZ 2 
Gt 

Pt - Pi 


2PL + Pr) 


V 


Gt 


■ Mi( p l + 2 p r) 


(A14) 
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APPENDIX B 

EFFECTIVE MOMENT OF INERTIA OF BULKHEADS 


The contribution of the cover sheets of the box beam to the effec- 
tive moment of inertia of the bulkheads can be approximated by an 
analysis of the plate -stringer combination shown in figure 17 . 

In the analysis the plate and stringer are assumed to be of dif- 
ferent materials but of constant cross-sectional dimensions. The con- 
tribution of the plate is expressed as an effective width which can be 
used to determine the area of an equivalent stringer having a total 
elongation under load the same as that of the stringer in the combined 
structure. The method of least work (page 156 , reference 6) is used to 
determine the state of stress in the plate and stringer. The stringer 
stresses can then be integrated over the length of the stringer to 
determine the total elongation. 


The stress distribution in the plate can be defined in terms of a 
stress function 0 as follows: 


dv = 


d 2 <J> 

dy 2 


a 2 $ 




xy 


d 2 0 
Sc 3y 


(Bl) 


The stresses given by equations (Bl) automatically satisfy equi- 
librium conditions; the differential equation which the stress function 
must satisfy to fulfill compatibility is 


dx 






Sx 2 dy 2 


= 0 


dy 


A solution of equation (B2) is given by 


(B2) 


<J> = 


n=l L - 


B n + c n (l + ^ 


n-Tty 


sin 


mtx 

l 


(B3) 
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In this expression Bn and Cn are arbitrary constants which are 
determined from the condition that the true stress distribution is that 
which makes the strain energy of the combined structure a minimum. If 
the plate width is assumed infinite, the stress function given by equa- 
tion (B3) satisfies the boundary conditions that all stresses vanish 
at y = 00 and that a x = 0 when x = 0 or Z. The stress function 
does not provide zero shear stress along the edges x = 0 and Zj 
therefore, in effect the plate has ribs along these edges. This viola- 
tion of boundary conditions is considered unimportant since, in the 
actual structure, these edges are restrained by the spar flanges. The 
strain energy in the plate is given by the expression 



- 2|ia x cry + 2(1 + 



dy dx 


(B4) 


Making substitutions from equations (Bl) and (B3) and integrating 
this expression gives 


it V B n E n c n c n 

P Z 2 n=l n \ 2G P 2G P + E P 


(B5) 


since 


c _ P 

G P - 2(1 + n) 


Consideration of equilibrium of plate and stringer gives the 
following expression for the load in the stringer 


p s = p o - (P 0 - p z)f - 2t JJ a x dy 

The strain energy in the stringer can he written as 


U 


s " 




dx 


(B 6) 


(bt) 


which, upon substitution from equation (B6) and integration, becomes 


U s = 


Z 

6ae c 


(p 0 2 + PoPz + p z 2 ) 


Jt 2 t 2 

AE 0 


n=l s n=l. 


(B8) 
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The total strain energy of the system is then 


U = Ur, + Uc 


(B9) 


If the total strain energy (equation (B9)) is minimized with respect 
to Bn and Cn, two equations are obtained which yield the following 
expressions for the constants 


B n = 


1 u 2 [po 

- (-i) n PzJ 

An 3 n 2 ( l + n)(3 _ n) 

^ E s _ kit 

E p An(l + p)(3 - n)_ 


C n ' M£) p 

- MH* 


The total extension of the stringer can now be determined as 

r l Pc 


u = 


AE S 


dx 


(BlOa) 


(BlOb) 


2AEs^ P ° + " 


2t_ 

AE S 


YZ fl - (-D^jBn ( B11 ) 

n=l 


The effective area of the equivalent stringer is defined as 

Ae = A + X,t 0 


(B 12 ) 


where 


te = *1? 


The elongation of the equivalent stringer is then 

r 1 P* 


u = 


0 ^e^s 


dx 


+ P 


l 


2E S vA + X,t e 


(B 13 ) 
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Setting equation (Bll) equal, to equation (B13) gives the following 
equation for the effective width 


where 



(bU) 


Note that the applied loads do not appear in the expression for 
effective width. 


A similar result can be obtained from an analysis in which the 
plate is assumed to be infinitely stiff in the transverse (y) direction. 
In that case, the differential equation of the plate is 


d 2 u + /G\ d^u 
bx 2 V E 'P by 2 


= 0 


(B15) 


The solution of equation (B15) which satisfies all the boundary condi- 
tions for a plate of width h is 


cosh 


u 




\l 


m* - y) 


G -n l 


n=l 


cosh, 


(IL 


sin 


nnx 

l 


(bi6) 


With the stresses defined in terms of u displacements as follows 

g - E ^U 

' E P s 


r r ^ 
r xy = Gp ^ 


(B17) 


and the strain energy of the plate given by 
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J£L 


dy dx 


The procedure previously followed can be used to obtain 



where 


(B18) 


(B19) 



Examination of equation (B19) reveals that the effect of finite 
plate width h is negligible whenever h I since the hyperbolic 
tangent terms very nearly approach a value of 1. In this case only a 
small difference is found between equations (Bl4) and (B19). 


The two expressions for effective width (equations (Bl4) and (B19)) 
have been evaluated for a range of values of A/lt e when \i = ^ 
and h = 00 . The results are plotted in figure 18. 


The results of the preceding analysis can be presented in a more 
convenient form for the analysis of a swept wing. The moment of inertia 
of the bulkhead alone can be related to the area of a substitute stringer 
as follows: 


I = |Ac2 


and the effective moment of inertia is 


I e = I + Ixtec" 


(B20) 


1(1 + v ) 


(B21) 
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where 



Figure 8 is a plot of v as a function of c 2 lt e /l, the results of equa- 
tion (Bl4) being used in the solution. 
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APPENDIX C 

NUMERICAL EXAMPLE 
Description of Specimen 


The application of the method to an actual structure is illustrated 
by an analysis of the untapered, 45° swept box beam of references 1 
and 2. Its construction details a.nd principal dimensions are shown in 
figure 9- The outer and carry-through sections are of doubly - 
symmetrical cross section and are divided into five and three bays, 
respectively, by internal bulkheads placed perpendicular to the spars. 

It is important to note that the actual structure contains a short 
bulkhead (bulkhead 7) within the triangular section and that its presence 
is ignored in the analysis because the method of analysis assumes a tri- 
angular section with no internal bulkheads. Figure 10 is a photograph 
of the specimen under test and illustrates the manner in which it was 
supported at the four corners of the carry-through section. In the 
analysis which follows these supports are assumed to be rigid. 

The dimensions of the three sections of the idealized structure 
are summarized in table II and illustrated in figures 19 and 20. The 
dimensions of the two bulkheads and that portion of the front spar 
bordering the triangular section are given in table III together with 
their calculated stiffness factors. The material of the specimen was 
24S-T3 aluminum alloy except for the steel bulkheads. These materials 
are assumed to possess the following elastic properties: 



E 

G 

Material 

(psi) 

(psi) 

24S -T3 

10.5 X 10 6 

4.0 x 10 6 

Steel 

29.0 X 10 6 

11.0 x 10 6 


Loading Conditions 


The box beam is analyzed under four different loading conditions, 
symmetrical and anti symmetrical bending and torsion corresponding to 
the test data of references 1 and 2. These loads are applied to the 
loading bulkhead at the tip of the outer section; the bending being 
produced by a vertical shearing force of 2,500 pounds; the torsion by a 
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pure couple of 43,420 inch-pounds. For convenience, these loadings 
are hereinafter referred to by number as follows: 


Loading 1 
Loading 2 
Loading 3 
Loading 4 


Symmetrical tip bending load of 2500 pounds 
Anti symmetrical tip bending load of 2500 pounds 
Symmetrical tip torque load of 43,420 inch-pounds 
Anti symmetrical tip torque load of 43,420 inch-pounds 


Analysis of the Outer Section 

Method of analysis .- The dimensions of the idealized outer section 
are given in table II and illustrated in figure 19 which also shows the 
notation that will be employed. The stress and distortion distributions 
for the idealized structure will be determined as well as* the stiffness 
factors required for the analysis of the complete structure. 

Since the outer section is a cantilever box beam on a flexible 
root , the stresses and distortions can be obtained from the superposition 
of the following solutions: 

(a) Outer section with a rigid root and a tip bending load of 

2500 pounds 

(b) Outer section with a rigid root and a tip torque load of 

43,420 inch-pounds 

(c) Outer section with the root warped 

(d) Outer section displaced as a rigid body 

A simplification of the analysis and the use of elementary theory 
in some instances are possible because of the constant doubly- 
symmetrical cross section. Further simplification is made by assuming 
that all bulkheads are rigid in their own plane, although this assump- 
tion leads to a slight violation of continuity because bulkhead 6 is 
assumed flexible in the analysis of the triangular section. Figure 21 
shows the two types of force groups which appear in this analysis. 

Shear flows 'are shown in addition to the concentrated flange forces. 

The sign convention used with the warping group is that of reference 7 . 

Bending of a cantilever .- The outer section is considered to be 
rigidly built in at the root and loaded by a central vertical shea.ring 
force of 2500 pounds at the tip as shown in figure 22(a). The internal 
forces and stresses are those of the F-group and can be expressed as 
follows: 


NACA TN 2232 4l 

F = £ = ^(L - x) = 178.57(89 - x) (Cl) 

CT B = a F = Me = Vc (L _ x) = 9T> oi(89 - x) (C2) 

q c B = q c F = ^ = 178.57 pounds per inch (C3) 

g 

t c B = r c F = 1- = 2289 psi (CM 


The beam bends without twisting and the deflection of the center line 
consists of two components, that due to flexure and that due to shear 
deformation of the webs. Thus, 

w 6 = w CT + w T (C5) 


where 


V a = g|Tj(3L - x)x 2 = 0.43994(267 - x)x 2 


(c6) 


w 


T 


V 

2Gct c 


x = 0.0005723x 


(C7) 


The distribution of stress and deflection is given in table IV(a). 


Torsion of a cantilever .- The outer section is considered to be 
rigidly built in at the root and loaded by a pure couple of 43,420 inch- 
pounds at the tip as shown in figure 22(b). The internal forces and 
stresses are those of the X-group plus the shear stresses required to 
equilibrate the torque and can be expressed as follows: 


T Xn T 
= — 


(C8) 


(Tt) 


c ,n 




c ,n 


T n 

2bc 


+ <1 


c ,n 


Tn X n T - X n _i T 
2hc + 2a n 


(C9) 
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/ . \ T T _ in , x ^n 

^ Tt )b,n " %,n " 2b c ^,11 " 2bc 


Y 1 Y J 

A n " A n-1 


2a 


n 


(CIO) 


The X-forces are statically indeterminate j however , they can be calcu- 
lated by means of the following recurrence relation (reference 7 ) which 
establishes continuity between bays: 


•^n^n-1 " (in + Pn+l) ^n + ^n+l^n+l “Jn^n + ^n+l^n+1 (Cll) 

where 


_a_ 1 /_b_ _c_\ 

^ 3AE + 5Ga\tb + t c / 

f = 1 /I. . _c_ 

"6AE BGa\tt t c 

1 - ^ A- - 

d oGbc \tb t c ) 

in which A is the effective area for warping stress. 

Once the X-forces are known, the twist of one bulkhead relative to 
the next can be determined as follows (reference 8): 


- *n-l) (C12) 


For the dimensions and stiffness parameters given in table II, the 
recurrence relation yields the following set of simultaneous equations 
for the X-forces: 


A 6 


n 


. fit') * 


( 

’ (§) n - 4J n(C 


where 


J = 


2b 2 c 2 

iL + JL 

*b tc 
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-3.57806 

0 . 575 H 

0 

0 

0 


r~ t 
X 2 


0 

0.57511 

- 3.57806 

0.57511 

0 

0 


x 3 T 


0 

0 

0 . 575 H 

-4.15317 

1.75718 

0 


x 4 T 

= 

0 

0 

0 

1.75718 

-4.72828 

1.75718 


y T 
x 5 


0 

0 

0 

0 

1.75718 

- 2 . 36414 


y T 

3 _ 


3296.924 


The solutions are 


x/ = 

-10.573 

x 3 T = 

-65.779 

xi T = 

-398.673 

x 5 - 

-920.751 

T 

X6 1 = 

-2078.917 


The distribution of stress and twist is given in table IV(b). 

Warping of the root . - The root of the outer section is warped the 
amount shown in figure 22(c). The internal forces and stresses axe 
those of the X -group which can be determined by the method of analysis 
used in the preceding section. The only change is in the boundary 
conditions with the warping of the root being related to the X-forces 
as follows: 

f 6 x 5 “ - p 6 x 6 W - I («*2 - h) !C13) 

The solutions of the system of equations are then 

x 2 w = -5612(02 - 03) 

X 3 W = -34915(02 - 03) 

X 4 W = -211615(02 - 03) 
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x 5 W = -488733(02 - 03) 
x 6 w = -11O3484(0 2 - 0 3 ) 

Table IV(c) summarizes the stress and twist distributions for 

(0 2 - 0 3 j = 1 X 10- 6 . 

Rigid-body displacements .- The outer section is given the rigid- 
body displacements shown in figures 22(d) and 22(e), These displace- 
ments do not affect the stress distribution but give the outer section 
a twist about its center line 




(Cl*) 


plus a deflection of that center line 

^ = | w 2 + |(02 + 03) x (C15) 

Superposition of solutions . - The complete stress and distortion 
distributions can be obtained by combining components from each of the 
preceding analyses. Since the basic expressions are the same for both 
symmetrical and anti symmetrical loads, the stresses and distortions at 
any cross section can be written as follows: 



Loadings 1 and 2 

Loadings 3 and 

crp 

cr B + a w 

+ a w 


cjB _ a W 

-a T - aW 

T f 

-T B - T W 

c 'c 

-T T . T W 

c c 

T r 

B W 

T c - T c 

-T T - T w 

c c 

T b 

_ W 
-T b 

T T T W 
_T b - b 

0 

eR + e w 

eP + e T + e w 

w 

w^ + w® 

wR 


(ci6) 


The sign convention for the stresses is that shown for the internal 
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forces of the outer section in figure 2. Positive deflections and twist 
move the front spar downward. Similar expressions can he written for 
the forces at the root cross section from which the stiffness factors 
given in table V can be determined by inspection. These are the stiff- 
ness factors required for the analysis of the complete structure. 


Analysis of the Carry -Through Section 

Method of analysis . - The dimensions of the idealized carry-through 
section are given in table II ard illustrated in figure 20, which also 
shows the notation employed. The stress distribution is determined 
along with the stiffness factors required for the analysis of the 
complete structure . 

Since the supports at the four corners of the carry-through section 
are assumed to be rigid (wq = W3 = 0 ), the stiffness factors associated 
with the deflection and twist of the ends in the plane of the end cross 
sections will not be required. The stresses can therefore be obtained 
from the superposition of the following solutions: 

(a) Carry-through section with the end cross sections rotated out 

of their original planes 

(b) Carry-through section with the end cross sections warped out 

of their original planes 

The doubly- symmetrical cross section permits considerable simpli- 
fication of the analysis since in such cases the end rotation is a 
result of the application of F-forces only and the warping is the result 
of the application of X-forces only (fig. 21 )j however, a few complica- 
tions are introduced because the splices in the center bay make it 
stiffer than the other two bays. The bulkheads are assumed rigid in 
their own plane despite the fact that in the analysis of the adjacent 
triangular section, bulkhead 8 is assumed to be flexible. 

Since the analysis depends upon whether the end distortions are 
symmetrical or antisymmetrical, the ratio R is introduced in order 
that general equations, applicable to both types of loading, can be 
written. Then: 


R = +1 when the loading is symmetrical 

R = -1 when the loading is antisymmetrical 

Rotation of the ends .- The ends of the carry- through section are 
rotated, symmetrically or antisymmetrically by the amount shown in 
figure 23(a). Equilibrium of internal forces requires that: 


F F F 10 “ F 9 f 9 " f 8 

^c^lO aio 

(C IT) 

The forces can now be related to the distortions by means of equations 
follows: 

(A13) and (Al4) as 

*9 - *8 - 2 (ISi)( f 8 + F 9 ) 

(C18) 

w 9 ■ W 10 ' ^9 a io ' 2 ( 2 Gt c ) 10 " 2 (6ei) 10 ( 2F 8 + F io) 

(C19) 

“8 - “9 ' <* 9*9 - 2 (2Gt c ) 9 * 2 (sei) 9 ( f 8 + 2f 9 ) 

(C20) 

and since 


„ 8 = w lx = 0 

(C21) 

F 10 F 11 w 10 ^10 011 
~ F 9 " f 8 ” w 9 ~ ~$9~ ~ "08 

(C22) 

08 = ^(^1 + ^3 sin A + 0 ^ cos A \ 

(C23) 


the following relationships can be obtained: 



HlOJ 

ja 10 + (1 - R)a^j + ^3 sin A + 02 cos a) 

1 - R 

+ (1 - R) (f-\ 

vc /10 vc /9 

/ca\ ( 2 +R)aio 2 /ca\ 

0 „ _ d - R ) 2 P 

a io + (1 ” R ) a io a 9 + 3 ^ 

G 

\EI / 10 3 Ul/9 


(C24) 
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U 7 


F 9 = F 8 


1 + ^2-(l - R) 
a 10 


F F 1 - R 

q c ,9 q c,10 8 a 10 + (1 - R)a 9 


(C25) 


(C26) 


The stress distribution corresponding to end rotations of the amount 

+ 1^2 sin A + ^2 cos A = 1 x 10”^ is given in table Vl(a), in which 



(C27) 

(C28) 


Warping of the ends . - The ends of the carry-through section are 
warped, symmetrically or antisymmetrically, by the amount shown in 
figure 23(b). Then, since 


„ ^10 Ai *10 "'ll 

- x 9 - x 8 - 'If ‘ 'W 

and 

^8 = §(% - ^2 sin A ' ^3 cos A ) 


(C29) 


(C30) 


the method of reference 7 can be used to obtain the following equations 
for the warping forces: 


P 9 X 8 - fc*X 9 + j^T = - ^2 sin A - cos A ) (C31) 

-f 9 X 8 + (p 9 + P 1Q - Rf 1Q )x 9 + (j 1Q - J 9 )t = 0 (C32) 


The torque in the carry- through section is statically indeterminate for 
antisymmetrical loadings but can be determined from the condition that 
one end does not twist relative to the other because of the rigid 
supports, that is: 


A (? 9 + A0]_q + A<9]_2. = 0 


(C33) 
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which can be written as follows by use of the equations of reference 8: 
4(1 - R)j 9 x 8 - 4(1 - R)(j 9 - j 10 )X 9 + 


Q 

The stiffness parameters f, j, p, and — have been previously 

Go 

defined in the section M Analysis of the Outer Section" of this appendix 
and their numerical values for the carry-through section are listed in 
table II. These values can be used to obtain the following sets of 
simultaneous equations and solutions: 



Symmetrical loads: 


3.349163 

-2.874949 0.0797693 

x 8 


1.75^1 -^3 sin A - 0 3 cos A j 

-2.874949 

3.727325 -.0388418 

x 9 

= 

0 

0 

0 .0336914 

T 


0 _ 


Xq = 1546392(^1 - ^3 sin A - 0^ cos a) 

X 9 = 1192759 (*l - ^3 Sin A - 03 cos a) 

T = 0 


Anti symmetrical loads: 


3.349163 -2.874949 0.0797693 


X 8 


1*75 (^1 - ^3 sin A - 0 3 cos A'j 

-2.874949 7.096181 -.0388418 


X9 

= 

0 

.638154 -.310734 .0336914 


T 


0 


Xq = 1509857(^1 - +3 sin A - 0 3 cos a) 

X 9 = 479367 (^1 - t 3 sin A - cos A j 

T = -24177260 (^ - ^2 sin A _ 0 3 cos A^ 

The stress distribution corresponding to end warpings of 

(^1 - ^3 sin A - 0 3 cos a) = 1 x 10"^ is given in table Vl(b), in which 
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k9 


W _ 1 / T n x n ~ x n-l \ 

c j n ^c,n\^ Bc 2a. n / 

W _ 1 / T n ^ ~ X n-l\ 

T b,n “ t b ^ n \2bc " 2a n / 


(C35) 
(C3 6) 
(C37) 


Superposition of solutions . - The complete stress distribution can 
be obtained by combining the components from each of the preceding 
analyses. The basic expressions are the same for all types of loading} 
thus, the stresses on any cross section can be written as follows: 

N 


(Jp = 
CJ R = 


a B + a W 


a B - a w 


T .p = -T 


T„ = T , 


B 

c 

B 


- T 


- T 


w 

c 

w 


T, = -T, 


w 


(C38) 


The sign convention for the stresses is that shown for the internal 
forces of the carry-through section in figure 2. Similar expressions 
can he written for the end forces from which the stiffness factors given 
in table VII can be determined by inspection. These are the stiffness 
factors required for the analysis of the complete structure. 


Stiffness Factors of the Beams 

The dimensions and stiffness factors of the two bulkheads and the 
portion of the front spar which bound the triangular section are given 
in table III. 

Since the nature of the normal -stress distribution in the 
triangular-section cover influences the effective moment of inertia of 
the front spar, two values are given. They were determined as follows: 

For bending stresses, I e was taken as one -half the moment of 
inertia of the doubly- symmetrical outer section. 


50 


NACA TN 2232 


For warping stresses, I e was determined from the effective 
warping area of the outer section; thus, 

I e = ^Ac 2 = |-(0.86 3)(7 ) 2 = 21.1 inches^ (C39) 

The effective moment of inertia of the bulkheads includes a con- 
tribution from the cover sheet of the box beam as determined from fig- 
ure 8. The equivalent thickness of the aluminum-alloy sheet acting 
with the steel bulkhead is given by the relationship 

he = t - Esheet = o.05o(^j = 0.0181 inches (C40) 

^bulkhead V 29 . 0 / > 


Triangular Cover Sheet 

The shear stiffness of the triangular cover sheet, which frequently 
appears in the general equations, has the following value: 

~ 2 ~ - x 10 6 )(0.050)(7) = 700,000 pounds (c4l) 

The shear stiffness per unit width also appea.rs and is 

= ^(700,000) = 23,333 pounds per inch (C42) 


The Systems of Equations and Their Solutions 

Sufficient data have now been obtained to permit evaluation of the 
coefficients of the matrix (table I); however, since the supports have 
been assumed rigid ' (wq = w^ = 0), it will be unnecessary to evaluate 
the coefficients of wq/b and w^/b. Furthermore, there are only seven 
unknown joint displacements, which require only seven equations. The 
logical equations to eliminate are equations ( 3 ) and ( 9 ) since they 
contain the support reactions Vq and V 3 , which are statically indeter- 
minate for the anti symmetrical loading conditions. The joint force V 2 
is zero because the loads are applied at the tips only.' 

A different set of equations will be required for each loading 
condition since the structure responds differently in each case. .Thus, 
the stiffness factors for the carry-through section are different for 
symmetrical and antisymmetrica.1 loads and the effective moment of inertia 
of the front spar and the loading terms are different for bending and 
torsion. The calculated coefficients are given in matrix form in 
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table VIII. Each block in the table has space for four values , one for 
each loading condition. If a single values is given, that coefficient 
is good for all four loading conditions; when two values are listed, the 
upper one is for loadings 1 and 2 whereas the lower one is for loadings 3 
and and if four values appear, they are for loadings 1, 2, 3, and 4 , 
respectively, when reading from top to bottom. 

Throughout the calculations, a large number of significant figures 
have been carried, more than are justified by the accuracy of the initial 
data; however, the extra figures were carried in order to obtain an 
accurate check on the numerical work when the calculated internal forces 
are substituted into the original nine equilibrium equations. The solu- 
tion of the equations and the calculation of internal forces frequently 
involve the differences of large numbers and the final results are apt 
to contain several significant figures less than the initial coefficients. 

The solution of the systems of equations yields the joint displace- 
ments given in table IX. Many methods are available for the solution 
of simultaneous linear equations; however, the method of reference 9 is 
recommended because of its many practical advantages. 

Calculation of Stresses 

The flange forces and shear flows in and around the triangular 
section can be obtained from the joint displacements by substituting 
them back into the force -displacement relationships, equations ( 13 ) 
to (35)- The shear stress is given by the shear flow divided by the 
sheet thickness; thus. 


T 



«* 3 ) 


The flange stress in the front spar and bulkheads is obtained from the 
flange force and the effective moment of inertia as follows: 


a 



(Ckb) 


The flange stresses in the outer and carry-through sections must be 
determined by summing up the various component stresses, equations (Cl6) 
and (C 38 )j since different effective areas are associated with the 
bending and warping stresses. The results of these calculations are 
given in table X. 

The stress distribution in the outer and carry-through sections can 
be obtained as described in the analysis of these sections, with joint 
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displacements substituted where necessary. The results are given in 
tables XI and XI I , respectively. 


Fisu*"® 12 illustrates the distribution of flange normal stresses 
and spar shear stresses as calculated in this example and compares them 
with the experimental data of references 1 and 2. The distribution of 
normal stress in the cover at three selected stations is similarly 
illustrated in figure 13. 


Calculation of Distortions 

The distortions of the outer section can be calculated by adding 
up a number of component distortions as provided in equations (Cl6). 

The deflections of the individual spars can be obtained from the 
deflection and twist of the box beam as follows: 

Wp = w + —0b 

(C45) 

W R = W - -0b 

The twist of the structure in a plane parallel to the flight path con- 
tains a component of the twist perpendicular to the spars and a 
component of the rate of change of deflection, but it is most easily 
calculated from the deflections of the individual spars. 

The calculated results are listed in table XIII and are graphically 
compared with the experimental data of references 1 and 2 in figure 11. 
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TABLE I.- COEFFICIENTS (a^) OF THE GENERAL MATRIX 
DEFINING THE JOINT DISPLACEMENTS 


Coefficient 

Formula 

a ll 

|T k l5 + k l6) " ( k ll + k l^) sec A + ( a + 

3)12 cos A 

a 21 

-(a + 3)^2 cos A 


a 3l 

( k 19 + ^0) ” (^7 + ^8) 


a Ll 

-3^2 cos A ” S]_2 cos A 


a.51 

Gtc 

~2b~ COS A 


a 6i 

(a. - 3) 12 cos A + 2 ^e 12 + sin A 


a 7l 

( k 23 + ^2l-)^13 sec A - ^13 cos 2A 


a 8l 

(k-Q + ^12) sec A + sin A 


a^l 

“J_f 13 ( k 23 + + ( k 27 + ^8)1 b sec A + 

^r( e i3 cos ^ ' *23) 


a 12 

(a + 3) 23 tan A 


a 22 

(a + 3)23 cot A 


a 32 

0 


a k2 

Gtc ~ . . 

2b~ & 12 cot A 


E52 

°23 Gtc . . 

6^ + 2lT cot A 


a 62 

(a - 3)23 tan A - ^|^-(ei2 + *23) 


a 72 

(a + 3)23 + B13 cos 2A esc A 


a 82 

-(a + 3)23 tan A - 


a 92 

-(a - 3)23 sec A - ^—^13 cos 2A - €23) 

esc A 
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TABLE I.- COEFFICIENTS (a^) OF THE GENERAL MATRIX 
DEFINING THE JOINT DISPLACEMENTS - Continued 


Coefficient 


Formula 


aU 

a 5L 

a 64 

a 74 

a 84 

a 9^ 


a 13 

|j k 15 ' k l 6 ) " ( k ll “ k l 2 )_ 

tan A - (a + P ) 23 tan A 

a 23 

-(a + 3 ) 13 cot A - (a + 3 ) 23 cot A 

a 33 

jj k 19 - ^ 0 ) - (^7 “ ^8)] 

sin A 

a 4 3 

0 



p 23 


a 53 

5 23 


a 63 

(a - p ) 23 tan A 


a 73 

( k 23 ” k 2l)^13 tan A - 00^3 

- (a + P ) 23 

a 83 

( k ll " k 12) tan ^ + ( a + 3)93 ^an A 

*93 

■Jjj.3 (^3 “ k 2l) + (^7 ■ ^ tan A - 


(a - p ) 13 cos A - (a - 3 ) 23 sec A 

& lk 

(a + P ) 12 sin A 


& 2k 

-(a + p).^ sin A - (a + p ) 13 esc A 

a .34 

0 



Gtc 
2b 

cos A cot A 


-3l2 sin A + 3 ^— 5i2 cos A cot A 

Gtc 
2b 


G"b c 

(a - P) 12 sin A 2~( e 12 + e 23) cos A 

G"fc c 

-Pl 3 sec A + 2^— cos 2A cot A 
Gtc 


cos A 


Gtc 


(a - 3) 13 - — (613 cos 2A - e 23 ) cot 


A 
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TABLE 

I.- COEFFICIENTS (ay) OF THE GENERAL MATRIX 

DEFINING THE JOINT DISPLACEMENTS - Continued 

Coefficient 

Formula 

a 15 

-(a + P ) 12 

a 25 

(a + 3) 12 

a 35 

0 

a L5 

1 Gtc ^ 

k 3 + a 12 + 2 b“ 6 12 

a 55 

. Gtc 

_k 7 + 2 b" 

a 65 

-^1 - 623 )^ 7 ^* "^ an A + ( a “ 3)l2 “ + € 23 ) ^an ^ 

a 75 

G"b c 

-k^ tan A + — - S-y cos 2 A sec A 

a 85 

1 GtC 4 . A 

-k^ - ^ tan ^ 

a 95 

^1 - € 23 )^!) sec A - -^p -( € ]_3 cos 2A - ^ 23 ) sec A 

a i 6 

( k 15 " k l 6 ) ' ( k ll " k 12 ) 

a 26 

(a + 3 )j _3 

a 36 

(k^ - k 2 o) cos A - (k 2 y - cos A 

a A 6 

, Gtc 0 

- k 3 - Jr Sl2 

a 56 

k 7 - 5i£ 

a 66 

(l - er^k^h tan A + ^^( c i2 + e 23) ban A 

a 76 

k 3 tan A + 8^3 ^ 1^23 - kg^j + ctQ _3 tan A - 

^9. 5-i 0 cos 2A sec A 
2 b x 3 

a 86 

k 3 + (&q - k-]_ 2 ) + tan A 

^6 

-(i - «23) k 7 b sec A _ |^l3( k 23 " + ( /k 27 " k 2 ^] b + 

(a - P )]_3 sin A + cos 2A - e 23 ) sec A 


NACA TN 2232 


57 


TABLE I.- COEFFICIENTS (a-y) OF THE GENERAL MATRIX 
DEFINING THE JOINT DISPLACEMENTS - Continued 


Coefficient 


Formula 


a l8 

a 28 

a 38 

a 48 

a 58 

a 68 

a 78 

a 88 


a, 


98 


a l7 

[J k l7 + k l8) " ( k l3 + k l^) 

a27 

0 

a 37 

[0*21 + *22) " 0*29 + k 3o) 

a iv7 

(a - P) 12 cot A 

a 57 

0 

a 67 

2(a - 3) 12 cot A 

a 77 

(kp^ + ^26)^1 3b sec A + (a 

a 87 

^k 13 + k^jb sec A 

897 

_ [fl3( k 25 + **26) + 0*29 + 1 
2(a - 3)^3 cos A 


b sec A 


0 

0 

0 


-(a - 3) 12 cot A 


' (a " P) 23 

■2(a - 3)ip cot A - 2(a - 3)93 tan A 


0 

0 


2(a - 3) 23 sec A 
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TABLE I.- COEFFICIENTS (a^) OF THE GENERAL MATRIX 
DEFINING THE JOINT DISPLACEMENTS - Concluded. 


Coefficient 


Formula 


a 19 

a 29 

a 39 

a L9 

a 59 

a 69 

a T9 

a 89 


a. 


99 


Q k 17 - k l8) " ( k l3 ~ k iJ] b sec A 
0 

[^1^21 - k2 2 ) - (kqy - k^o) b 
0 

(a ' P) 23 623 
2(a - 3)23 tan A 

( k 25 " k 26) 5 13 t sec A " ( a " P)i3 
^23 - k-^j b sec A 

- ^13(^5 ' ^6) + (^9 “ k 3o)_ b2sec A 

2(a - P )]_3 cos A - 2(a - P)23 sec A 


a 10 


- (k]_ + kl^M 

^0 


(ki + kjM 

a 30 


■ i(vi + V 2 + 

a 4o 

k^M + k2T 

a 50 


-kgv - k^T 

a 6o 


* (*5 

a 70 

(k-|M - k2T ) 

a 8o 

1 

t-^M - k^T 

a 90 


"c V 3 + ( k l0 


b tan A 


b sec A 


TABLE II,- DIMENSIONS AND STIFFNESS PARAMETERS OF IDEALIZED 


STRUCTURE USED IN THE NUMERICAL EXAMPLE 


Item 

Outer section 

(fig. 19) 

Triangular 

section 

Carry -through section 
(fig. 20) 

Bays 2 , 3 and 4 

Bays 5 and 6 

Bays 9 and 11 

Bay 10 

Dimensions 

B, in. 

30 

30 

30 



b ' } in . 




lip I 

12. 1 

c, in. 

7 

7 

7 

7 

* 7 

tb, in. 

0.050 

0.050 

0.050 

0.050 

0.0885 

t c , in. 

0.078 

O.O78 

0.078 

0.078 

0.078 

a, in. 

22 

11 

30 

9.60 

9.76 

I f in^ 

90.2 

90.2 


. 122.58 

135.15 

**A, sq in. 

O.863 

O.863 


0.96k 

1.229 


Stiffness parameters 


***p, in. /lb 

1.78903 X 10" 6 

2.36414 X 10 ' "6 


3.349163 X 10 -6 

2.062590 X 10 ■ 6 

f, in. /lb 

0.57511 x 10" 6 

1.75718 X 10 


2.874949 x 10' 

1.684428 x 10~ 6 

S, per lb 

O.O75931 x 10 " 6 

0.075931 X 10 


0.0797693 x 10 

0.0409275 x 10~ 6 

a/GJ } per in. -lb 

0.0430112 x 10‘ 6 

0.0215056 X 10 


0.0128741 X 10 

0.0079432 x 10 -6 


Includes an allowance for splice plates. 

^Effective area for warping stress ^A = Af + ^ct c + ^btft + \ } A< 


•***, 


These stiffness parameters are defined in appendix C. 


VJI 
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TABLE III.- DIMENSIONS AND STIFFNESS FACTORS 
OF THE BEAMS 


Beam 

Front spar 1-2 

Bulkhead 2-3 

Bulkhead 1-3 

Material 

24S-T3 aluminum alloy 

Steel 

Steel 

Stress 

Bending 

Warping 



Dimensions 

l, in. 

30 

30 

29.25 

41.34 

c, in. 

7 

7 

7 

7 

t, in. 

0.078 

0.078 

0.125 

0.125 

I , in. 4 



6.15 

8.01 

t e , in. 



0.0181 

0.0181 

c 2 lt e /l 



4.21 

4.57 

V 

] 



0.60 

O.65 

Ie, in. 4 

45.1 

21.1 

9.84 

13.21 

Stiffness factors 

per lb 

0.575170 -xio 6 

0.743315 X 10 6 

0. 338627 XlO 6 

0.332182 X 10 6 

a, lb 

3.993615 XlO 6 

2. 400311 X 10 6 

4 . 346797 x 10 6 

4 . 427535 x 10 6 

P, lb 

0.516384 XlO 6 

-0.290336 X 10 6 

-l. 559397 xio 6 

-1. 780133 X 10 6 

S, in. 

11.14498 

8.623871 

4.295415 

4.51467 

€ 

0.257002 

0 . 425063 

0.706297 

0.781584 
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TABLE IV.- STRESS AND DISTORTION DISTRIBUTION IN OUTER SECTION 

(a) When root is rigidly built -in and tip bending load of 
2500 pounds is applied 


n 

X 

(in. ) 

f b 

(lb) 

B 

a 

(psi) 

t c B 

(psi) 

v a 

(in. ) 

w T 

(in. ) 

B 

w 

(in. ) 

1 

88 

176 

97 


0.6098 

0.0504 

0.6602 

2 

66 

4107 

2231 

2289 

.3852 

• 0378 

.4230 

3 

44 

8036 

4365 

2289 

.1899 

.0252 

.2151 

4 

22 

11964 

6500 

2289 

.0522 

.0126 

.0648 

5 

11 

13929 

7567 

2289 

.0136 

.0063 

.0199 

6 

0 

15893 

8634 

2289 

0 

0 

0 


(b) When root is rigidly built-in and tip torque of 
43,420 inch pounds is applied 


n 

a T 

(psi) 

-c T 

(psi) 

T 

T b 

(psi) 

Ta/GJ 
( radians ) 

AD* 

(radians) 

a T 

(radians) 

1 

0 





-o.oo684o 






2 

-13 

-1328 

-2063 

-0.001868 

0.000003 

-.004975 

3 

-76 

-1341 

-2043 

-.001868 

.000017 

-.003124 

4 

-462 

-1422 

-1916 

-.001868 

.000101 

-.001357 

5 

-1067 

-1630 

-1593 

-.000934 

.000159 

-.000582 

6 

-2409 

-2000 

-1015 

-.000934 

.000352 

0 


(c) When root is warped by amount 


(02 “ 0 3 ) = 1 X 1°' 6 


radians 


n 

aW 

(psi) 

T W — 

(psi) 

■, b u 

(psi) 

A0W 

' 1 ans ) 

e w 

(radians) 

1 

0 

-0.0016 




2 

-0.0065 

0.0025 

0.001705 

.333450 

3 

-.0404 

-.0085 

• 0133 

.008900 

. 324550 

4 

-.2452 

-.0515 

.0803 

.053668 

.270882 

5 

-.5663 

-.1615 

.2519 

.084167 

.186715 

6 

-1.2787 

-.3583 

.5589 

.186715 

0 
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TABLE V.- OUTER -SECTION STIFFNESS FACTORS 


Stiffness 

factor 

Loadings 1 and 2 

Loadings 3 and 4 

kjM, lb 

15892.86 

0 

kgT, lb 

0 

-2078.91 

k^, lb/radian 

-110348k 

-1103484 

k^M, lb 

15892.86 

0 

ktjV, lb/in. 

-178.5714 

0 

kgT, lb/in. 

0 

156.0246 

kj, lb/in. /radian 

27943 

27943 

k gV, lb/in. 

0 

0 

k^T, lb/in. 

0 

50.7371 

ki 0 V, lb/in. 

-178.5714 

0 




TABLE VI.- STRESS DISTRIBUTION IN CARRY-THROUGH SECTION 


(a) When ends are rotated by the amount 

sin A + cos Aj = 1 X 10“^ radians 


(b) When ends are warped by the amount 

sin A - 02 cos A) = 1 X 10 " 6 radians 


Stress 

Symmetrical 

Anti symmetrical 

B 

a Q > Psi 

1.3100 

0.7303 

B 

o , psi 

9- 

1.3100 

0.2460 

„ B 

a , psx 

9+ 

1.1882 

0.2231 

B 

T 0 , Psi 
c,9 

0 

1.6177 

B 

r , psi 

c,10 

0 

1.6177 


Stress 

Symmetrical 

An t i s ymme t r i c al 

w 

a 8 , psi 

l. 6 o 4 i 

1.5662 

W 

a , psi 

9 - 

1.2373 

0.4972 

w 

a , psi 

9 + 

0.9705 

0.3900 

W 

T n ' pSl 
c,9 

-0.2361 

- 1.2140 

_ w 

T , psi 

c,10 

0 

0.1038 

W 

T > Psi 

b ,9 

0.3683 

0.2530 

w , 

t , psi 

b,10 

0 

-1.0185 


ON 
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TABLE VII.- CARRY-THROUGH SECTION STIFFNESS FACTORS 


Stiffness factor 

Loadings 1 and 3 

Loadings 2 and 4 

kn, lb/radian 

3277276 

1827064 

ki2, lb/radian 

-15^6392 

-1509857 

lb/radian 

3277276 

1827064 

k , lb/radian 

1546392 

1509857 

k-^, lb/in. /radian 

0 

126178 

kgo> lb/in. /radian 

l84l8 

94691 

kg-^, lb/in. /radian 

0 

0 

kg^, lb/in. /radian 

-18418 

-12651 

kg -j, lb/in. /radian 

0 

-126178 

kgg, lb/in. /radian 

18418 

94691 
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TABLE VIII.- SYSTEM OF EQUATIONS FOR EACH LOADING CONDITION 


Loading 1 
Loading 2 
Loading 3 
Loading 4 


Symmetrical tip bending load of 2500 lb 
Anti symmetrical tip bending load of 2500 lb 
Symmetrical tip torque load of ^ 3,^20 in. -lb 
Anti symmetrical tip torque load of 43 , ^-20 in. -lb 


Loading 1 
Loading 2 
Loading 3 
Loading 4 



Loadings 1 to 4 


7.562904 

7.459566 

5.865823 

5.762485 

2.787400 

-5.880185 

- 5.807114 

-5.880185 

- 5.807114 

3.189065 

1.491984 

- 4.510000 

- 2.109975 

-3.092785 

- 3.019714 

-3.092785 

- 3.019714 

0 

-3.189065 

-1.491984 

2.787400 

- 5.434802 

- 6. 933040 
- 5-235959 

4.510000 

2.109975 

2.647402 

0 

-0.549024 

0.063013 

0.260050 

0.201223 

0 

- 0.181257 

0.347587 

5.357149 

3.705018 

- 1.363534 

- 1.304707 

-3.477231 

-2.690647 

-0.016499 

1.035295 

0.363038 

0.016499 

- 0.004610 

0.004610 

-1.375000 

2.935556 

2.462581 

5.231885 

5.114242 

5.906194 

1.981974 

1.342586 

2.556713 

1.652486 

0.920518 

1.038161 

-18.766850 

-17.193682 

-0.117596 

-0.080792 

-0.117596 

-0.080792 

2.787400 

-7.131781 

-7.157806 

-7.131781 

-7.157806 

2.517477 

- 1.103484 

5.614173 

5.588148 

5-614173 

5.588148 

0 

2.942805 

0.943574 

2.942805 

0.943574 

- 3.487400 

7.611069 

6.124321 

7.611069 

6.124321 

-0.494977 

- 1.803484 

6.627153 

5 . 1 4 0 405 
6.627153 

5 . 140405 

0 


^3 


0i 


02 


0 3 


wg 

T 


-0.031786 

0 


0.031786 

0 


0.015893 

- 0.0020789 


0 

- 0.0000507 


- 0.005357 

0.0057558 


0.015893 

0.0020789 


0.015893 

0.0020789 
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TABLE IX.- JOINT DISPLACEMENTS FOR EACH LOADING CONDITION 


Displacement 

Loading 

1 

Loading 

2 

Loading 

3 

Loading 

4 

radians 
^ 2 , radians 
radians 

0.0003681 

0.0029693 

0.0012872 

0.0019304 

0.0040942 

0.0024238 

- 0.0002970 

-0.0013392 

-O.OOO 6535 

-O.OOO 5616 

- 0.0015269 

-0.0008414 

01 _, radians 
02 , radians 
03 y radians 

0.0005019 

0.0055292 

0.0038611 

0.0005217 

0.0066615 

0.0050147 

-0.0003813 

- 0.0012307 

0.0001279 

-O.OOO 3817 

-0.0014186 

- 0.0000604 

wp/b, radians 
W 2 /b^ radians 
w^/b, radians 

0 

0.0025716 

0 

0 

0.0037002 

0 

0 

-0.0011405 

0 

0 

- 0.0013282 

0 
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TABLE X.- CALCULATED LOADS AND STRESSES AROUND THE TRIANGULAR SECTION OF 
THE SWEPT BOX BEAM FOR EACH LOADING . CONDITION 


Loading 1 
Loading 2 
Loading 3 
Loading b 


Symmetrical tip bending load of 2500 lb. 

Anti symmetrical tip bending load of 2500 lb. 
Symmetrical tip torque load of 43,^20 in. -lb. 
Antisymmetrica.1 tip torque load of 43,420 in. -lb. 



Flange 




Web or 

cover 







p n 



a n 




In 



T n 






(lb) 



. (psi) 



(lb/in. ) 



(psi) 



Loading 

Loading 

Loading 

Loading 

n 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 

1 

8110 

8148 

-1005 

-1006 

44o6 

4426 

-1167 

-1168 

87.7 

86.9 

78.5 

78.5 

1124 

1114 

1007 

1007 

2 

14052 

14076 

-580 

-580 

7634 

7647 

-673 

-673 

-132.0 

-132.6 

118.1 

118.1 

-1692 

-1699 

1514 

1514 

3 

17734 

17710 

580 

580 

a l0767 

a l0739 

672 

672 

225.2 

224.6 

118.1 

118.1 

2887 

2879 

1514 

1514 

A 

4689 

4656 

-1911 

-1911 

11675 

11593 

-4758 

-4758 

-46.6 

-46.0 

88.7 

88.7 

-932 

-920 

1774 

1774 

5 

2342 

2349 

-1365 

-1365 





-110.4 

-110.7 

64.3 

64.3 

-2208 

-2214 

1287 

1287 

6 

18197 

18164 

-2307 

-2306 

a i050i 

a lo466 

a -996 

a -995 

-60.3 

592.0 

1.4 

-144.1 

-773 

7589 

18 

-1848 

7 

8077 

8110 

-2076 

-2076 

&2 

a 37 

a -756 

a -757 

60.3 

42.1 

-1.4 

-1.0 

1205 

843 

-28 

-19 

8 

9224 

9231 

1761 

1761 

17107 

17120 

3266 

3266 

-60.3 

-1222.5 

1.4 

158.5 

-773 

-15673 

18 

2032 

9 

-5737 

-5764 

710 

711 

-10640 

-10690 

1317 

1319 

292.6 

311-5 

26.1 

25.7 

2341 

2492 

209 

206 

10 

14052 

14076 

-580 

-580 

a 6501 

a 6528 

-672 

-672 

219.6 

219.4 

-39.5 

-39.5 

1757 

1755 

-316 

-316 


£L 

Indicates a stress composed of a bending and warping component. 
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TABLE XI.- CALCULATED STRESS DISTRIBUTION IN THE 
OUTER SECTION OF THE SWEPT BOX BEAM FOR 
EACH LOADING CONDITION 


FLANGES 


n 

(T-p 

(psi) 

a R 

(psi) 

Loading 

Loading 

1 

2 

3 

4 

1 

2 

3 

4 

6 

6501 

6528 

-672 

-672 

10767 

10740 

672 

672 

5 

66 22 

6634 

-298 

-298 

8512 

8500 

298 

298 

4 

6091 

6096 

-129 

-129 

6909 

6904 

129 

129 

3 

4298 

4298 

-21 

-21 

4432 

4432 

21 

21 

2 

2220 

2220 

-4 

-4 

2242 

2242 

4 

4 


WEBS 


n 

T f 

(psi) 

T r 

(psi) 

Loading 

Loading 

1 

2 

3 

4 

1 

2 

3 

4 

6 

-1691 

-1699 

1513 

1513 

2887 

2879 

1513 

1513 

5 

-2020 

-2023 

i4n 

1411 

2558 

2555 

1411 

1411 

4 

-2203 

-2204 

1352 

1352 

2375 

2374 

1352 

1352 

3 

-2275 

-2275 

1329 

1329 

2303 

2303 

1329 

1329 

2 

-2286 

-2286 

1327 

1327 

2292 

2292 

1327 

1327 


COVERS 


n 

T b 

(psi) 

Loading 

1 

2 

3 

iT 

6 

-932 

-920 

1774 

1774 

5 

-420 

-413 

1935 

1935 

4 

-134 

-132 

2025 

2025 

3 

-22 

-22 

2061 

2061 

2 

-4 

-4 

2066 

2066 
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TABLE XII.- CALCULATED STRESS DISTRIBUTION IN THE 
CARRY -THROUGH SECTION OF THE SWEPT BOX BEAM 
FOR EACH LOADING CONDITION 


FLANGES 


n 

CTp 

(psi) 

a R 

(psi) 

Loading 

Loading 

1 

2 

3 

4 

1 

2 

3 

4 

8 

2 

36 

-756 

-757 

10500 

10466 

-996 

-995 

9 - 

1202 

112 

-784 

-257 

9300 

3424 

-968 

-333 

9 + 

1587 

305 

-722 

-238 

7939 

2903 

-866 

-298 


WEBS 



T f 

(psi) 

T r 

(psi) 

n 

Loading 

Loading 


1 

2 

3 

4 

1 

2 

3 

4 

9 

10 

-773 

0 

7589 

11977 

18 

0 

1 1 

H H 
vog 

00 00 

-773 

0 

-15673 

-11284 

18 

0 

2032 

1932 


COVERS 



T b 

(psi) 

n 

Loading 


1 

2 

3 

4 

9 

1205 

843 

-28 

-19 

10 

0 

-3391 

0 

77 
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TABLE XIII.- CALCULATED DEFLECTIONS AND ROTATIONS OF THE OUTER SECTION OF 
THE SWEPT BOX BEAM FOR EACH LOADING CONDITION 


Loading 1 
Loading 2 
Loading 3 
Loading 4 


Symmetrical tip bending load of 2500 lb 
Anti symmetrical tip bending load of 2500 lb 
Symmetrical tip torque load of 43,4-20 in. -lb 
Anti symmetrical tip torque load of 43,420 in. -lb 


Deflection 

(in.) 


Rotation 0 

(radians) 


Front spar 


Rear spar 


Perpend i cular 
to spars 


Parallel to 
flight direction 


Loading 


Loading 


Loading 


Loading 


3 


0 

11 

22 

44 

66 

88 


0.077 
• 153 
.252 
.507 
.818 
1.159 


0.110 

.200 

• 311 

• 591 
.928 

1.293 


-0.034 

-.053 

-.072 

-.112 

-.152 

-.192 


•o.o4o 

-.061 

-.082 

-.126 

-.170 

-.215 


0 

.067 

.161 

.414 

.725 

1.065 


0 

.080 

.187 

.464 

.800 

1.166 


.006 

.014 

.029 

.045 

.061 


.005 

.010 

.021 

.033 

.045 


0.00257 

.00288 

.00302 

.00311 

.00313 

.00313 


0.00370 

.00401 

.00415 

.00423 

.00425 

.00425 


-0.00114 

-.00198 

-.00287 

-.00471 

-.OO657 

-.00844 


-0.00133 

-.00216 

-.00305 

-.00489 

-.OO676 

-.00862 


.00556 

.00751 

.00860 


•00557 

.00752 

.00862 


.00206 

.00336 

.00468 


-.00206 

-.00336 

-.00468 


^Positive deflection is downward. 

^Positive rotation is front spar downward. 
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Figure 2.- Exploded view of idealized structure showing internal forces 
on each component around the triangular section. 
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Figure 4.- Positive directions of joint displacements. 
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P L- 


■ P R 
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(a) Loads. 




(b) Distortions. 



(c) Differential element. 


Figure 5-- Loading and distortions of beam. 
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(a) Assumed stress distributions in actual structure. 



(b) Equivalent flange forces in idealized structure. 

Figure 6.- Assumed distributions of stresses and forces in actual and 

idealized structures. 
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(a) Actual cross section. 
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'S3 
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(b) Equivalent cover. 


b 



(c) Idealized cross section. 


Figure 7*- Relationship between actual and idealized cross sections. 
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I 

Figure 8.- Chart for determining the effective moment of inertia of 

bulkheads. I e = l(l + v) . 
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Sheet splice 
Flange splice 



Section A-A 


Figure 9*- Details of sweptback box beam used in tests and analyzed in 

illustrative example. 



Figure 10.- Anti symmetrical bending test setup of sweptback box beam. 


s NACA r 

L-5758^ 



Deflection, in. Deflection, in. 
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(a) Symmetrical tip bending load. 

\ 


2 Experimental (ref. 1,2) 

Theory without shear lag 

Theory with shear lag 



(b) Antisymmetrical tip bending load. 

Figure 11.- Experimental and theoretical deflections and rotations of 

sweptback box beam. 


Deflection, in. t Deflection, in. 
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(c) Symmetrical tip torque load. 


§ Experimental (ref. 1,2) 

Theory without sheor log 

Theory with shear lag 



(d) Anti symmetrical tip torque load. 


Figure 11.- Concluded 
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(a) Symmetrical tip bending load. 



(b) Anti symmetrical tip bending load. 


Figure 12.- Experimental and theoretical spanwise distribution of flange 

normal stresses and spar shear stresses of sweptback box beam. 

* 
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(c) Symmetrical tip torque load. 


0 Experimental (ref. 1,2) 
— Elementary theory 

Theory without shear lag 

— Theory with shear lag 



(d) Antisymmetrical tip torque load. 

■ 

Figure 12.- Concluded. 
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(a) Symmetrical tip bending load. 


0 Experimental (ref. !, 2) 



(b) Anti symmetrical tip bending load. 

Figure 13*- Experimental and theoretical chordwise distribution of normal 

stress of sweptback box beam. 


Stress, ksi Stress, ksi 
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(c) Symmetrical tip torque load. 


° Experimental (ref. 1,2) 



43.42fm.-kips 


43.42 in.- kips 



(d) Anti symmetrical tip torque load. 


♦ 


Figure 13*- Concluded. 
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* 



(a) Stress distribution 
in actual structure. 



(b) Flange forces in 
idealized structure. 


Figure lk . - Shear-lag stresses and forces in actual and idealized 

structures. 



Figure 15- - Type of idealized structure required to include shear-lag 

effects in the analysis. 
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(b) Symmetrical torsion. 

Figure l6.- Calculated effects of bulkhead flexibility on the stresses 

and distortions of a swept box beam. 
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t Figure 17 • - Plate -stringer combination analyzed to determine the effective 

width of cover acting with bulkhead. 



Figure 18.- Effective width of plate acting with stringer. (See 

appendix B) . 
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Figure 19.- Idealized outer section used in illustrative example. 



Figure 20.- Idealized carry-through section used in illustrative example. 
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Figure 21.- Force groups of outer or carry-through sections used in 

illustrative example. 
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(a) 


(b) 


R=+l for symmetrical loads 
R=-l for antisymmetrical loads 



Figure 23.- Types of individual analyses of carry-through section used 

in illustrative example. 
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